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This report describes work performed by Software and 
Engineering Associates, Inc. (SEA) under NASA contract 
NAS8-34974 "A Shock Wave Capability for the Improved 
Two-Dimensional Kinetics (TDK) Computer Prograr". The TDK 
computer program is a primary tool in applying the JANNAF liquid 
rocket engine performance prediction procedures. The purpose of 
this contract has been to improve the TDK computer program so 
that it can be applied to rocket engine designs of advanced 
type. In particular, future orbit transfer vehicles (OTV) will 
require rocket engines that operate at high expansion ratio, 
i.e. in excess of 200:1. Because only a limited length is 
available in the space shuttle bay, it is possible that OTV 
nozzles will be designed with both relatively short length and 
high expansion ratio. In this case, a shock wave may be present 
in the flow. An objective of the present study has been to 
modify the TDK computer program to include the simulation of 
shock waves in the supersonic nozzle flow field. The shocks 
induced by the wall contour can produce strong perturbations of 
the flow, affecting downstream conditions which need to be 
considered for thrust chamber performance calculations. Project 
manager for this project has been Mr. Gary R. Nickerson. Dr. 
Lanh Dang has performed the computer programming and has 
provided engineering support. The project was very much aided 
by the helpful support of the contract monitor, Mr. Klaus W. 
Gross, and by Mr. A. Kresbach. 
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1 . INTRODUCTION 


The Two-Dimensional Kinetics (TDK) 1 computer program is a 
primary tool in applying the JANNAF liquid rocket thrust chamber 
performance prediction methodology. 2 This computer program and 
the per .■ rmanee prediction methodology were originally developed 
under the auspices of the performance Standardization 
Subcommittee (PSS) of the JANNAF. The goal of the PSS is the 
development of a methodology that includes all aspects of rocket 
engine performance from analytical calculation to test 
measurements* that is physically accurate and consistent, and 
that serves a 3 an industry and government reference. 

Recent interest in rocket engines that operate at high 
expansion ratio, such as most OTV design, has required an 
extension of the analytical methods used by the TDK computer 
program. The primary objective of the study reported here has 
been to equip the TDK computer program with the capability of 
simulating a shock wave induced by curvature of the nozzle wall 
contour. This capability includes the effects of chemical 
reaction kinetics. For example when Hydrogen and Oxygen are 
used as propelants, the dissociation of water behind an oblique 
shock front is simulated. 


In order 
streamline to 
to revise the 
has been made 
l .e . mass a 


to accomodate variation in mixture ratio from 
streamline within the flow, it has been necessary 
transonic analysis method used by TDK. No attempt 
to treat mixture ratio variation along streamline, 
ddition or diffusion effects. No attempt has been 


Ref 


Nickerson, G. R., Coats, D 
•'The Two-Dimensional Kinet 
Program,” Engineering and 
sytems , Inc., December 197 
N A F 9-12652, NASA JSC. 


. E. , and Bartz , J . L. 
ic (TDK) Reference Computer 
Programming Manual, Ultra- 
3, prepared for Contract Nc. 


Ref. 2. "JANNAF Rocket Engine Performance Prediction and Calcu 
1 at i on Manual", CPIA 246, April 1 975. 
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made to treat the interaction o r shockw;ves with slip 
discontinuities. Instead, the approach has been to allow a 
continuous variation in mixture ratio to be specified transverse 
to the flow, and to trace by construction methods the 3 hock path 
and flow interaction. 


The work carried out in the study is summarized below: 

a) In the TDK program, the existing clip line 
method has been converted to a streamline 
method to avoid unrealisti- shock interaction 
with the slip line. Each streamline is 
assigned a specific mixture ratio. Completely 
coupled two-dimensional flow with finite rate 
chemical kinetics is computed. 

b) The existing TDK transonic modei^ has been 
modified to analyze flow with variable 
mixture ratio. 

c) The method of character i sti c (MOC) concept 
has been modified to utilize the varying 
properties of the streamlines. 

d) Various techniques for constructing the 
supersonic flow field with the MOC concept 
have been investigated, and several shock 
construction methods have been examined. The 
most efficient ones were selected for 
incorporation intr the TDK program. 

e) A set of supersonic flow control subroutines 
were written to construct the supersonic flow 
field with shock waves included. 

f) The TDK program with shocks was checked out 

for flows with gas properties along 

streamlines that are either: 

1 ) constant , 

2) chemically frozen (i.e. fixed com- 
position, but properties varying with 
temper atur e ) , or 

3) governed by finite rate kinetics. 

g) Thu TDK program, coupled with a boundary 
layer module, was adjusted such that the 
effect of a shock wave undergoing regular 
reflection at the nozzle wall was accounted 
for . 

h) T- iniques were introduced for reducing 
c . puter time . 
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METHOD OF APPROACH 


The TDK computer program described in Reference 1 is the 
JANNAF approved computer program used for the purpose of 
calculating two-dimensional rocket nozzle flows. The program is 
highly accurate and contains many useful features. However, no 
provisions have been made for including the effects on nozzle- 
performance caused by the occurrance of strong shock waves. 
Previously, it had been assumed that a nozzle wall which induces 
a strong shock in the flow should be considered a poor design 
and, thus, not a design of interest. More recently, nozzles 
containing shock waves have become important. There are several 
reasons for this. The foremost reason is that the engines 

required for Orbit Transfer Vehicle (OTV) propulsion feature 
unusually high expansion ratios (A/A*>200). In an effort to 
shorten the length of these nozzles, designs are being proposed 
that induce shocks in the flow. If the shock strength is 
everywhere small, then the flow streamlines remains isentropie 
and the TDK program will accurately compute the nozzle 
performance. The basic reason for this is that the shock 
strength, 6, defined as: 

6 ^ (P 2 - 

effects a change in entropy in accordance with the following 
well-known expression: 

AS/R - y + 1 6* * higher order terms. 

1 2 y 2 

Thus, weak shocks are essentially isentropie. 

According to Hoffman (Reference 3. Volume II, page ) , the 
method used by TDK will yield "reasonable*' results even for 
relatively strong shocks, i.e. shocks such that 6c 1 *. 

p e l 71 — Zuorow, M.J., and Hoffman, J.D., Gas Dynamics , Volume I 
and II, John Wiley & Sons, Inc., New York, 1976. 


Nevertheless , since OTV Mission requirements are highly 
sensitive to th- engines delivered specific Impulse, It Is 
crucial that TDK he able to accurately compute the effects of 
sh< ek waves on performance, Nozzle shock waves are Inherently 
two-dimensional In their character, so that It la necessary that 
a two-dimensional analysis be used to compute the, effects, 
-mce TOK IS the approved JANNAF program for calculating 
effects of two-dimensional phenomenon on nozzle performance (see 
Reference «5 . It Is appropriate that It be extended to Include 
snook waves. The method of approach used for achieving this 
objective is presented below. 

n is required that TDK be able to accurately assess the 
effects on performance of shock waves that are Induced by . 
nozzle wall. This accuracy requirement was an over: 1 ng 

consideration when selecting the approach to be used. Thus, e 
numerical solution selected is based on the Methcd of 

Characteristics (HOC), since It represents the post accurate 
method available. Shock tracing (as opposed to shock capturing) 
was also .selected for reasons of accuracy. 

Although accurate, the MDC with shock tracing Is not 
suitable for flows that feature many shocks, since the logic to 
treat the reflection, Intersection, and coalescence of m.ny 
shocks is intractable. It was assumed, however, that only a ew 
strong shocks are present, and that th.se originate from he 
wall. Multiple regular reflections from the flow axis and 
nozzl e wall (see Figure la) were treated. Mach reflections from 
the flow axis or the nozzle wall (see Figure lb) were not 
treated. It was assumed that Mach shocks, which Include 


Ref .4 


"JANNAF" Rocket Engine * 

evaluation Manual," CFIA Publics. 246, Aprii, 
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t il> urc 1 : Reflection of iin oblique klimt wuve from « »•«■ *ui Kevlar 

reflection ib i Mach reflection 


region of subsonic flow, will not effect nozzle performance. 
T:.at is, if a Mach shock occurs, it is assumed that its’ domain 
of dependance does not include the nozzle wall. 

in the present study, it is also assumed that a shock 
originates by the coelescence of Right Running Characteristics 
( RRCs ) due to flow angle gradient change, or turning, at the 
nozzle wall. The case where the shock starts attached to the 
wall, due to an abrupt discontinuity in wall slope, is not 
treated at this time. Only the more difficult induced shock 

case is treated. 

Since the shock wave is very thin, it is assumed that the 
chemical composition across the shock will be frozen. On the 
downstream side of the shock, the reacting flow will relax 
towards a new equilibrium state at increased pressure. If the 
flow is near equilibrium, this process can occur extremely fast, 
causing numerical difficulties that can be severe. This problem 
has been solved by utilizing a fully implicit numerical 

integration method. 


* 


The version cl' TDK tru t is ueacr i bed in Reference b was 
taker, a; the starting point for the .resent fograrn . All of the 
capabilities described in Reference 5 have been included in 
the newer program. 


Ref. 5 Nickerson, G.R. and Dang, L.D., "Improved 

Two-Dimensional Kinetics (TDK) Computer Program", SEA 
Inc., Final Report, SN56, NASS-BS^Ob, October 1983. 
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3 . TECHNICAL DISCUSSION 

The technical effort is discussed under the f 
riven below; 

Transonic Model with Mixture Ratio Variation 
Selections of Shock Methods, 

Shock Construction Methods, and 
Techniques for Reducing Computer Time 


3.1 1 ran: cni e Mode l with M 1 / 1 r ■ I • : r. • 


The transonic model used by "C K is described in References 1 
and 6. it is a f irs t-order perturbation analysis that in 
based on the method of Sauer . Toe method divides the flow into 
regions of constant mixture ratio that are separated by 
sliplines. Each rf.’ion contain.., a specified fraction of the 
total mans flow. This method, which is called the striated flow 
option, has been moaified so that flown with mixture ratio 
variation can be analyzed. This variable mixture ratio option 
is described below. 

When the variable mixtu; e ratio option is used, there are no 
sliplines in the flow. Instead, the flow mixture ratio ,r, will 
vary from the axis (¥-0) to the wall (v-l) as specifiec by an 
input table of r versus S'. The streamline function .S', 
represents the mass flew between the streamline and the axis, 
divided by the total nozzle ma3S flow. An ODK calculation is 
done for each entry in the above table. Value. along the 
initial data line for the MO C are obtained by interpolation ir, 
the ODK results sing pressure and radial coordinate position 
,Y, as independent variables. The transonic analysis is used 
to provide a table of S' versus Y. The method used is described 
below. 

The ODK program constructs tables of flow properties 
(p, V, T, and c ) as a function of pressure. These tabl es span 
the nozzle throat region. An average expansion coefficient is 
computed using these tables as 

Y - In (P./P )/ln (p./p ) 

*’ i * i 


Ref. 6 Nickerson, G. R., "Striated Flow in a Converging- 
Diverging Nozzle", Dynamic Science Report CS-2/71-1, 
prepared for NASA JSC, February 1971. 

Ref. 7 Sauer, R. "General C har a ct er i s t i cs of the Flow Through 
Nozzles at Near Critical Speeds", NACA Tech. Note No. 
111)7 (19H7). 

8 

& V.i r- ' ' ' 


l 


"T" 


where the subscripts 1 and £ refer to tae first and last table 
entries, respectively. Values of Y are found from Y and from 
input to the transonic analysis 

Y n _ 1 )/2 n • 1.2, ...N 

and also 


5 - * - ¥ n _, 

*n n n-i 

Using these values, the transonic analysis calculates 

Y - 0, Y , Y , . . .Y N - Y wall * 

0 12 

The above T„ values represent the radial location at which the 
input mixture ratios 

r , r , r , . . . .r N 
0 1 * 


are located. In this way the input table of r versus * is 
converted to a table of r versus Y. 


the transonic analysis 
3 0 • B 1 ’ C 1 ’ C 2_ 


Next , 

coefficients (B~ , B. . C. , C„ ). 


is used to compute the 
These are used to compute 

P(X,Y) and 6 (X,Y) in the transonic region (see Reference 1, 
pp.2-20) at points n-0,1...N. Using eacn of these N + 1 values 

of P as an independent variable, the corresponding values for 

the 


p , V , T , and c i 


are obtained by linear interpolation from 


corresponding table that was computed by ODK. These tables are 


^ !*-. -v^v ' ;■ 


9 


then used to linearly interpolate for P,p,V,S,T, and c . at each 
MOC initial line point using ^ as the independent variable. 

The program will not function properly if the spacing in the 
mixture ratio table, r R , is too large. The required spacing 
depends on the chemoial system. As a rule each entry must 
differ no more than A or 51 from its adjacent values, depending 
on the stoichiometry of the system. 


The average engine mixture ratio, r aye , is also calculated: 


wall 


wall 


f (JL.) dm / [ ( — ) dm 

1 r+1 } r +1 


where 


r is the mixture ratio at posi 4 ion Y, and 


,v iLlLllil v d v 


3.2 Selection of Shock Methods 


Several techniques were investigated for constructing the 
supersonic flow field using the Method of Characteristics (KOC) 
concept. The MOC construction methods considered were: 

1) Combinations of RRC , LRC , and Streamline 
tracing , 

2) Streamline Normal method, 

Q 

3) Hartree scheme . 

It was determined that streamline tracing must be used 
because of the numerically stiff behavior of the finite rate 
chemistry. Streamlines must be traced and a f ully-ifoplicit 
integration method must be used in order to accurately compute 
the reacting flow chemistry while maintaining a reasonable 
integration step size. 

Since the shock waves to be analyzed come from the nozzle 
wall, it is necessary to trace Right Running Characteristics 
( RRC * s ) . A shock wave is inititiated when adjacent RRC's 
cross-over. Since computing the effect of these shock waves on 
nozzle performance is a primary objective of the study, a RRC 
construction method was considered to be necessary. Since the 
Streamline Normal Method and the Hartree scheme use 
interpolation formulas to locate characteristic lines, those 
methods were considered to be less accurate than RRC tracing. 

It was concluded that Method 1, (combinations of RRC, LRC, 
and Streamline tracing), must he used if the effect cf an 
oblique shock on engine performance is to be accurately 
estimated . 

Two shock construction methods were investigated: 

a) Shock capturing, and 

b) Shock tracing. 

It was concluded that Method a, above, must be used, again 
for reasons of numerical accuracy. 

Ref .8 Fox, L., Numerical Solution of Ordinary and Partial 

Differential Equations, Pergamon Press ,1962 . 


1 ] 


3 . 3 Shock Construction Methods 


In order to construct supersonic nozzle flow fields 
containing, shock waves, a set of logic control subroutines and a 
set of point calculation procedure subroutines were developed. 
The flow field point calculation subroutines ares 


SH0CK 

INPTR 

SHCKR 

SHCKA 

SHCKA1 

SHCKL 

SHCKW 

SHCKW1 


Basic running interior point 
Right running interior point 
Right running shock point 
Axis shock reflection point 
Off-axis shock reflection point 
Left running shock point 
Wall shock reflection point 
Off-wall shock reflection point 


The flow field construction procedure is outlined below. 


The mesh construction begins at the intersection of the 
initial data line and the nozzle wall and i-.RC’s are constructed. 
After the nozzle wall is reached, succesi ”2 RRC’s may cross. If 
so, a shock is inserted into the flow field at the cross-over 
point. Next, LRC's are constructed starting at the wall, and 
tne region up to the cross-over point is filled in. The LRC 
construction then continues with the first point on each LRC 
being a right running shock point. When the axis is reached, 

left running shock. The program 
construction scheme and the shock is 


traced, but 


shock i 

s r 

ef lected 

as 

reverts 

t 0 

a RRC 

cc 

d until 

it 

reaches 

the 

and is 

val 

culated 

as 

ds the 

fl 

ow axis 

• 

pie refl 

ect 

ions are 

al 


An example calculation is shown in Figure 2, which presents 
the mesh construction for a source flow with 1 I/ 1 * degree cone 
half angle that enters a right regular cylinder. A Right 
Running Shock (RRS) wave occurs at the intersection .f the cone 
and cylinder. In Figure 2 the shock is attached to the wall. 
If the induced shock opt is used (rather thar the attached 
shock option) a shock w A *x be detected when the first RRC from 
the cylinder crosses the RRC that is immediately upstream. The 
control subroutine detects this crossing and inserts a RRS 
point. The RRS is traced to the axis, reflected, and then 
traced as a Left Running Shock (LRS) to the wall, reflected, 
etc. Four axis reflections and three wall reflections are shown 
in Figure 2. 

Strictly speaking, regular reflection from the axis 
cannot occur in axiSymmetric flow. This is because the term 
(sin e)/r in the MOC relations becomes infinite behind the first 
ihock front; i.e., r is zero, but 0 is not zero. Thus a Mach 
shock will occur. However, as a practical matter, unless the 
shock is strong, the Mach disk is very small and can be ignored. 
In the computer program the Mach disk is ignored by reflecting 
the shock based on assuming that the incident shock angle at the 
axis is the same as at the previously calculated off axis shock 
point. 

The flow axis can be thought of as being a right cylinder 
with zero radius. If the radius of this cylinder were non-zero, 
regular reflection would be possible only if the radius exceeded 
some minimum value. At each RRS point this condition is 
checked, i.e. the Mach number is calculated behind a reflected 
shock resulting in zero flow angle. The value must be greater 
than one for the calculations to continue. 





th* ee 


Figure 3 shows Mach number vs distance for 
streamlines: the wall streamline, the median streamline, and the 
axis. It can be seen that the shock is much stronger at axis 
than at the wall. The strongest shock occurs at the first axis 
reflection. The Mach number behind the reflection is 1.58. A 
very rapid steepening of the shock occurs near the axis, such 
that in actuality, there is a Mach disk at the axis with 
subsonic flow downstream. This small region is ignored by the 
computer program. 

Figure 4 shows pressure versus distance for the same 
streamlines as shown in Figure 3. It can be seen from these two 
figures that there is a weakening of the shock front as the flow 
progresses down the cylinder. Eventually, the flow in the duct 
will be at constant pressure. 


The shock wave is constructed using the point calculations 
listed at the beginning of this section. Details of two of the 
subroutines are presented below. The basic shock wave 
calculation, subroutine SH0CK, is used to compute properties 
across an oblique shock front that is inclined at an angle, 6, 
with respect to the approach streamline. Subroutine SH0CK is 
described in Section 3.3.1. The Right Running shock point 
calculation, subroutine SHCKR , is described Section 3.3.2. It 
is used for shock initiation, as well as for RRS point 
cal culation . 
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Figure 3: Mach Number versus Distance for 3 Streamlines 




3.3.1 Subroutine SHOCK 


Given the shock angle, £• , and the properties ft point 
1, this subroutine finds all values at point 2 behind the 
oblique shock front, set figures 5a and 5b, below. 


(Step 1) 

The quantities u, v, and M, are calculated as 

u 1 = sin e 
Vj = Vj cos 6 

M 1 = v l / (> 1 RT 1 ) 1/2 . 

The Mach number, M 1 , is used to obtain a first estimate for 
P ^ 2 ^ es follows (equation 128 of NACA 1135): 

p( 2 )= P l{ 2Y l M l sin2s - ( \ - 1) }/ (Y 2 4 1) . 

(Step 2) 


Using the relations, 
: 1 U 1 = ° 2 U 2 


F + C, U 2 = P_ + p„u 2 


1“1 


2 2 


iVi 


°2 U 2 V 2 


5 


] •m'' 

* tr; +h. 


Mass Conservation (1) 

Normal Momentum (2) 

Tangential Momentum 
Energy Conservation (3) 


P/P • RT , R = '-\ /M,, Gas law 

w 

and 


M) 


h * h(T) 


Gas Properties (5) 


IS 





The following are calculated: 
from (1) and (2) 


'2 


(P 1‘ P 2 + ^u^/PjUj 


f t om ' 3 ) 


h 2 = \ “l + h l " 5 “2 


from ( r ) using subroutine TCALC 


T 2 * T 2 <V 


from (] ) 


C 2 * »l V U 2 


from (4) 


F^ i+1) = P 0 R T 2 


Iterate on P^ until 

|p< i+1) - e' 1 ’ 


/ P 


(i) 


< € . 


(Step 3) 

Values at point 2 are calculated as follows: 
6 * arc tan (v^/u^) - arc tan (7,/u,) 


1' u r 


e, + 6 


V 2 * ^2 + ^ 2 ] 1/2 


All other properties are obtained from Subroutine GPF. 

If M 2 is returned from GP* such that 

M. < 1 + c , nominally e = .2, 

2 m m 

an error message is printed and the run terminated. 
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3,3.2 Subroutine SHCKR 

This subroutine calculates a pair of right-running 
shock points, 3a and 3b, using the known points la, lb, 4, 
5, 6, and 7 shown in the figure below. 



The procedure used consists of 8 steps, as follows: 

1) 'Jse point 5 if available, otherwise 
calculate from la and 4, using INPTRS. 

2) Using 4 and la, including 8 1 , locate point 3 

* 

on line 4-5 . 

3) interpolate along 4-5 for all properties at 
3a. 

4) Use shock subroutine to calculate properties 
at 3b from 3a. 


p c x the first shock point that is computed, 
where us arc sin 1/M evaluated at the average 
back point and point la. 


set 8 * 
of the folded 



5) Interpolate for po.nt 2 between points f. and 
7, Set poin - . 6 os lb for first estimate. 

6) Relocate point 3 by intersection of RRC line 
2-3b and streamline 4-3a. 

7) Calculate g from slop of line 3-3. 

P) Return to step 3 and iterate until 
convergence . 

Details of this procedure are given below. 

Step 1) 

Point 5 is calculated (unless available) from points la and 
4 using subroutine JNPTR just before entering this 

subroutine . 

Points la, lb, 4, 5, 6 and 7 are known points and 

points 3a and 3b are the unknown points to be computed. 
Initial values (i=0) at point 3a for A, B, 0 , Y , and 6 are 
taken as the values at point 

Double subscripts imply averaged values, e.g. 



Step 2 ) 

Using points 4 and la, point 3 is located on the line 4-5. 


(i) r l ~ x l tan (P l +6 la ) " r 4 + X 4 tan 6 < r 
x 3 = tan e 45 - tan ( B i +e ia> ~~ 

'3 1 ’* ‘A + (x 3 1> ' *«> ta " 6 « 

Step 3) 

Properties at point 3a are determined by interpolation and 
integration along the line 4-5, as follows; 

C * U 3 -x 4 ) / U 5 -x 4 ) 


e (i+i) . e 4 + c <e 5 -e 4 ) 


, 0 + 1 ) . 
3a 


,(i+l) 
3 a 


P 4 + C <P 5 -P 4 ) 


, ( i+1 ) 

3a 


&1 


i) 


43a 


exp 


•( 


h — )* (x<l + 1 >-x 4 ; 

cos e / 43a \ 3a 4/ 


„ ( i + 1 ) _ T 

*2a < 


, ( i + 1 ) 
3a 


ny 


(i) 


i y x 


'4 3a 


exp 



— s_y 11 (m 

cos e ) 43a \ 


x (i+1) ->- I 
x 3a x 4l 


Next the gas velocity is calculated as 


N 


o+i) ln ) 

' in RT 1 ^! 


2N 


(i+1) 


v (i + l) = |v 2 4 - - — 

V 3s T'l N <1+1 )_ ! 0 ( 


N (1+1) -l 


1 > 


,(i + D 1 

3*' 


N 


(i+1) 


1/2 
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The species concentrations at point 3a , cj * * are 

calculated using the Species integration Subroutine, SINT, 

,(i+l) 

3a 




. , ( i + 1 ) R (i+i) 

and the gas properties at point 3a, a 3 „ • D 3a 

_(i+l) „(i+l) o (l+1) Y (i+1) , are calculated using the 

G 3a ' H 3 s ' 3a ' T 3a 
GPP subroutine. 


Step 4) 

Calculate all properties at point 3b from 3a using the 
shock wave subroutine, SHOCK. 


Step 5) 

Values at point 2 are determined by interpolation 
between points f end 7. Points 6 is initially assumed to be 
point lb and point 7 is the next point on the back LHC data 

line . 


The radial location of point 2 is calculated from 


(i+ 1 ) 


r< i+1, -[^ r 6 ♦ * r> ] 

frk^ tan i [6 u» : 4>'. «“>] 

l 7 


a 3b] 


If r ^ > r , point 7 becomes point 6 and t.ie next point on 
the back data line is point 7. The above calculation is 
repeated until r g < r^ +1 ^< * 7 * 

The axial location, flow angles, and Mach angle of 
point 2 are calculated from 


C 


2 


r < i+1) . r 

r~ - r. 


6 


(i+ 1 ) 

2 


+ c 2 lx 7-*6> 
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e ,W) 

* e 6 + C 2 <W 

^ 1+1) 

“ “6 + C J < W 


Step 6 ) 

Relocate point 3 by the intersection of RRC line 2-3b and 
streamline 4-5 


x (i+l) 

3 


-r 2 + x 2 tan ^ ( e 2 ~ « 2 + e 3b - a 3b ) + r^ - x^tane^ 
tan \ (e 2 - a 2 + e 3b - a 3b ) - tan e 45 


r (i+ 1 ) 
r 3 


r 4 + 


(x^ 1+1) - x 4 ) tan e 45 


Step 7) 

The revised shock angle is calculated from 

= ‘ e 13a * arc tan 

B 3 = 2 3 " 6 1 

where | 6 3 | > <* 3 
Step 8 ) 

Return to Step 3. 

The above integration equations are iterated (i^l,...) 
until r 3 , X3 , 63 » P 3 , P 3 # T 3 , and V ? converge to the 
required accuracy. 


(r l “ r 3 i+1) )/(Xj - x^ i+1) ) 
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Techniques for Reducing Computer tIiup 


3 .« 


work was carried out on the ODK and TDK nodules for the 
pose of reducing computer run time. Three aeparate tasks 
we;: carry out 1, single precision in ODK .’reus 
PPOOI -on. 2) modification of the matrix inversion method, and 

i ;:::?;v he chemicai — - — 1: 

ST01CC . Each task is discussed below. 


3 .«. 1 


ODK Preclsl on 


It was found that double precision is not needed In t e 
linear e,uatlon solver. DESK. when using DEC 32 bit word 
computers. The program was run with and without double precision 
for test case with high expansion ratio. The results showed 
no appreciable differences. The change from double precision to 
single precision In ODK saved about 30J of execution time. 


TDK Matrix Inversion Procedure 


NESK t JT llnear aqUat ‘° n 5 ° 1Ver USea by TDK is subroutine 

vector 5 T' !! S ° 1Va thS linCar e ’ Uati0 " -T. the 

... ' ' Chemical species mass fractions. In an 

thatT 1 1° reaU ° a C ° mPUter the Pnocedure was changed so 

that * was found and saved for use on sucessive Iterations. 
Two methods were investigated. 


In the first method a" 1 
elimination and saved for successive 


is computed by Gaussian 
iteration. This method has 



the advantage of being conceptually clear. However, it pays o: f 
only if there are at least some minimum number of iterations, 
depending on the size of tne linear system. The approximate 

number of operations for the AX»b method and for the a” 1 method 

> I 2 

are, respectively, kn and 5n + (k-1)n , where n is the size 

IF? 

linear system (i.e., the number of species) and k is the number 
of iterations. Thus, for a system of 10 species, at least 4 
iterations must be made before execution time is reduced. 

A second method was investigated in which the sequence of 
row exchanges and row operations in the Gaussian eliminations i 

were saved the first time, and the same sequence of operation ' 

were applied to b on successive iterations. This method is f 

slightly more complex than the previous one, but has the 
advantage that execution time is always equal to or less than • 

that of the straight Gaussian elimination at 'every iteration. 

Morever, additional savings in execution time can be achieved by 
using scaled Gaussian elimination instead of unsealed Gaussian 
elimination. ' l 

i : 

The following table shows the timing of the different ( 

methods for our test case. \ 

. t I 

* i 

ME THOD TIMING (fee.) * SAVING , 


Present 

Met hod 

164 


Method 

1 

162 

1 .2* 

Method 

(unsealed 

2 

Gaussian ) 

158 

3.7* 

Method 

2 

152 

7.3? 


(scaled Gaussian) 



The drawback of the new methods is that they have not yet 
taken into account the slight dependence f A on the step size, 
which is variable. Although this adjustment can be done, the 
added complexity and risk involved outweighs the small saving in 
execution time. Thus, it was decided to leave the linear 
equation solver in TDK the way it is. 


3.^.3 Modifications to Subroutine ST0ICC 

In order to reduce the number of calls to subroutine 
ST0ICC, the stoichiometric coefficients are obtained at the 

beginning of the run and then stored in various arrays. With 
this change, the subroutines DERI V , SDERIV, EE and EF2D do not 
need to call ST0ICC repeatedly and, therefore, the execution 
time is reduced. In our test cases the saving was about 5?. 
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. Conclusions 






4 





.1 

• r 



The TDK computer program that is described is Reference 5 
has been modified to Incorporate the work items a) through h) 
listed in Section 1 . All of the capabilities described ir 
Reference 5 have been retained. 

The TDK program has been given the capability to detect an 
oblique shock wave induced by curvature of the nozzle wall. The 
shock is detected by a crossing of right running chracteristics . 
The shock is traced through the flow field, including multiple 
regular reflections from the nozzle axis and wall. Only one 
such shock wave is allowed. Mach shocks are not allowed. 


The chemistry of the flow with a shock can be either a) 
constant properties, or b) frozen chemical composition or o) 
finite rate kinetics. Equilibrium chemistry with a shock wave 
is not allowed. TDK has not been provided with a exit plane 
option for flows with shock waves. 


Straited flow (i.e. dividing the flow into regions of 
constant mixture ratio separated by sliplines) is not allowed 
ith a shock wave. This i3 because shock-sl i pi i ne interaction is 
not provided in the program. The program can, however, be run 
with a continuous variation in mixture ratio from streamline to 

streamline. The transonic analysis has been modified to provide 
this capability. 


Modifications were made to the computer program to reduce 
computer run time and also to improve accuracy. Only a modest 
(on the order of 5 %) improvement in run time was achieved. 
Program accuracy was increased by several means. Complete 
energy conservation was achieved along streamlines in the HOC 
module by computing the gas velocity directly from the energy 
equation. The Bernoulli equation was used for a first estimate. 
Conservation of mass flow was forced by adjusting the flow 
pressure across each complete cnar acter i sti c surface. The gas 
density was then adjusted using the gas law. These measures 
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served to conserve energy flux and momentum flux across the 
flow. Accuracy was also improved by distributing points along 
the MOC initial data line using a sinusoidal distribution 
function. This metho:> provided an Improved MOC grid spacing. 

A revised User's Manual for TDK (i.e. Section 6 of Reference 
1) has been prepared. This document has been attached as 
Appendix A. The revised computer program and manual render 
obsolete the earlier versio-as-presented in Reference 5. 
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6.0 


PROGRAM USER’S MANUAL. 


ascription of the Computer Program Input. 

The TDK computer program consists of five modules, ODE, ODK, TRANS, MOC, and 
LM. All of these modules are required to perform a complete two dimensional 
ion-equilibrium nozzle performance calculation with a boundary layer. Various 
options exist in the program, however, which exercise the above modules alone, in 
part, or in combination. 

Data is read by the program sequentially in the order required for the execu- 
tion of the modules. This order is as follows: 

Tnermodynamic data, 

Data common to the modules, 

ODE module inputs, 

ODK module inputs, 

TRANS module inputs, 

MOC module inputs, and 
BLM module inputs. 

A more detailed description of these input data sets is presented in Table 6-1. 
The documentation in which each of the data sets is completely described is also 
.ndicated in Table 6-1. 

Of the data sets listed in Table 6-1, only the SDATA data set is required for 
, ,ery computer run. Input of the other data sets is required only if the options 
-„ney contain are to be used. For example, input of the thermodynamic data is not 
-equired if an existing thermodynamic data file is to be used. 
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Table 6-1. Input Data Sat Description 


Card Input Section No., 


Description 


THERMO 


thermodynamic data cards, see Tables 6-2 and 6-3 


END 

LOW t cphs 
• 

6.1.1, 

extension of thermodynamic data to temperatures 
below 300°K, See Table 6-5 

• 

r iD LOW T CPHS 

title 

6.2, 

one or more title cards 

DATA 

6.3, .. 

data directive card 

SDAIA 

SEND 

6.3, 

data namelist for module selection, and geometry 

REACTANTS 

6.4, 

reactants directive card 

• 

6.4.1, 

reactants cards., see Tables 2-6 and 2-7 

• 

6*4.1 f 

blank card required to end reactants cards 

OMIT 

6.4.2, 

cards to omit species 

INSERT 

6.4.2, 

cards to insert species 

kamelists 

6.4.3, 

0DE directive card 

SODE 

6.4.3, 

0DE namelist 

"PECIES 

• 

6.5.1, 

species cards for 0DK 

• 

r _ACTIONS 

6.5.2, 

reaction cards for 0DK 


LAST REAX 
l:r£RTS 


6. 5. 2. 6, inerts cards for 0DK. (to inert those species that 
are not named in reactions) 
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Card Input 


Section No.. 


Description 


THIRD BODY 6.5. 2. 7* reaction rate ratios for third bodies 


LAST CARD 


$ODK 

6.5.3 

0DK module namelist 

SEND 

STRANS 

6.6 

TRANS module namelist 

SEND 

SMOC 

6.7 

^5C module namelist 

SEND 

SBLM 

6.8 

ELM module namelist 

SEND 




Tbble 6-1 can be used as a guide when preparing input for given problem. It 
lists the data sets in the order in which they must appear in the data deck, and 
also shows the special cards which must appear in each set (first card, last card, 
etc.) if the program is to function properly. The table is basically self-explana- 
tory when used together with the detailed input descriptions which follow. 

Certain special options to the computer program are described separately in 
Section 6.9. 

An input data card listing for a sample case is presented in Section 7, follo- 
wed by the corresponding computer output. In preparing input to the computer prog- 
ram it is useful to review this input card listing. 

Successive cases can be run using the computer program but complete data should 
be input for each case. 
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6.1 


Thermodynamic Data . 


Ordinarily, a thermodynamic date, file is available for use with the program, 
and is assigned to logicial unit 25. The input described here can be used to gene- 
rate a thermodynamic data file if one is not available. 

This data set is identical to the THERMO DATA described in Appendix D of NASA 
’d _273 (i.e. Reference 9). 

Using this data set, thermodynamic data curve fit coefficients may be read from 
;ards. The curve fit coefficients are generated by the PAC computer program desc- 
ribed in. NASA TN D-4097 (i.e. Reference 22). 

The thermodynamic data (i.e. C® , etc.) are expressed as functions of tempera- 

ture using 5 least squares curve fit coefficient (a^) and two integration 
constants (a* -,) as frllows: 


+ a 7 

(a^y and a^) are specified for 
upper respectively. For the data 
interval is 300° to 1000°K and the 


C° 

P T 

T“ 


a 1 + a^T + a-T 2 ♦ a^ ♦ a-T 4 


H° 


a-T a,? 2 a.T^ a-T* 


h a-T a,T^ a.ir a c ur a* 

ic- a i ♦ -f ♦ ♦ 4" ♦ 4" * r 


'°T • a.T 5 a-T 4 

TT- *1 1«T ♦ V ♦ 4- ♦ -f- *4- 


For each species, two sets of coefficients 
two adjacent temperature intervals* lower and 
available in Reference 9 the lower temperature 
upper temperature interval is 1000°K to 5000°K. 


Ref. 22. McBride, B.J. , and Gordon, S. , "Fortran IV Program for Calculation of 
Thermodynamic Data", NASA TN-D-4097, Aug. 1967. 


The Input format required for this thermodynamic data ia defined in Table 6-2. 
ZT. “"k , f0r ^ SP9ClSS W ’ H ’ ”2°' M 2* °' OH. and 0, are Hated in Table 

"UDDlled^uith^h mody,iamic coefficients for many chemical apeciee are 

-applied with the computer program. A liat of theae apeclea la preaented in able 

0-4 # 


Data Tape Generation and Usage: 

A computer fun uaing thermodynamic data card input will generate a data tape on 

-og cal unit JANAF. This tape may then be saved and used at a later time. The 

program writes the THEHMO data card Images on unit JANAF as read but with two minor 

exceptions. The THERMO code card and the card numbers in card column 80 are 
omitted#- — — 


If ttarmodynamic data oarda are not input, the program assumes the theno- 
ds™mlCfdata ia on logical unit JANAF. Logical unit JANAF la currently assigned a 


TABLE 6-2 


FORMAT FOR THERMODYNAMIC DATA CARDS 


Card 

order 

Contents 

Format 

Card 

column 

1 

THERMO 

3A4 

1 to 6 

2 

Temperature ranges for 2 sets of co- . . 
efficients: lowest t, common T, and . 
highest T 

3F10.3 

1 to 30 

3 

Species name 

3A4 

1 to 12 


Date 

2A3 

19 to 24 


Atomic symbols and formula 

4(A2,F3.0) 

25 to 44 


Phase of species (S,L, or G for solid, 
liquid, or gas, respectively 

Ai 

45 


Temperature range 

2F10.3 

46 to 65 


Integer 1 ... 

115 

80 

4 

Coefficients a.' (i ■ 1 to S) in equations 
(6-1) to (6-3) 1 

5(E15.8) 

1 to 75 


(for upper temperature interval) 
Integer 2 

15 

60 

5 

Coefficients in equations (6-1) to 6-3) 
(ai, ai for upper temperature interval 
ana aj , ag, and a^ for lower) 

5 (El 5. 8) 

1 iu 75 


Integer 3 

IS 

80 

6 

Coefficients in equations (6-1) to (t-3) 
(a^,a 6 ,a 7 for lower temperature in- 

4 (El 5. 8) 

1 to 60 


Integer 4 

120 

80 

(a) 

r epeat cards numbered 1 to 4 in cc 80 
for each species 



(Final 

card) 

END (Indicates end of thermodynamic 
date) 

3A4 

1 to 3 


8 Gaseous species and condensed species with only one condensed phase can be 
in any order. However, the sets for two of more condensed phases of the same 
specie; must be adjacent. If there are more than two condensed phases of a 
species their sets must be either in Increasing or decreasing order according 
to their temperature intervals. 
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tABLE 6-3. THERMO DATA CARDS F.OR. AN 0 2 /H 2 PROPELLANT 
(Species AR, H, H 2< H 2 0, N 2 , O, OH, end 0 2 ) 


THERMO 

300,000 
A* 

0.250000006 

*0,t45375O26 

0. 

H 

0,250000006 

0.254716276 

0* 

M2 

0.310019016 

-«*B773S042E 

0.5S210391E 

H20 

0,271676336 

-0,299058266 

3.2963740AE 

Ht 

0.2896319AE 

•O.90S8M64E 

•0*632175596 

0 

0*254205966 
0*292308036 
*0* 160284326 
OM 

0«29106427E 

0,393538186 

0.187139726 

02 

0.362195386 
-0,120198286 
-0,676351376 
END 


1000,000 5000,000 
L S/66AR 100 


01 

03 


000 000 06 3l00,000 


01 0 

08-0 

0 

01 0 

03-0 

-08-0 

01 0 

05 0 
>08 0 

01 0 

03 0 
-09-0 

01*0 

05 0 
•08 0 

01 0 

04 0 
-09-0 

01 0 

04 0 
.08 0 


436600066 01 0*250000006 
*0 *.745374986 
J 9/6SH 100 000. 000 

460117636 00 0.250000006 
0*254716276 

J 3/61H 20 00 00 

SU19464E-03 0.526442106* 
19629421E 01 0.308744516 
181227396-1 1-0.968904746 
J 3'6lH 20 . 100 000 

294513746-02-0 *802243746* 

663056716' 01 0,407012756 
807021036-12*0,302797226 
J 9/6SW 20 00 00 

151S4866E-02-0 .572352776* 
616151486 01 0*367482616 
225772536-12-0.106115886 
J 6/6?0 10O 000 000 

27550619E-04-0, 31 0280336* 
492030806 01 0*294642876 
369069646-12 0*291476446 
J 3/660 1H 100 000 

959316506-03-0,194417026- 
544234456 01 0*383759436 
225710946-12 0.364126235 
J 9/650 20 00 00 

736182646-03-0*196522286- 
361509606 01 0*362559856 
215559936-1 Uo, 1^4752266 


* °* 
01 0 , 


5000.000 

0. 

0. 

03 0,436600066 01 

00 300.000 5000*000 

0, 0,. 

01 0, 0* 

08-0,460117626 00 

06 300,000 5000*000 

•07-0,349099736-10 0*369483486-14 

01 0,267652006-02-0,580991 626-05 

03-0,229970566 01 
0G A 300.000 5000,000 

>06 0,162266826-09-9,484721456-1* 
01-0,110844496-02 0,415211806-08 
05-0,322700466 05 
06 300,000 5000*000 

06 0,998073936-1 0-0, 652235556-14 
01*0 *120815006*02 0.232401026-05 

04 0,235804246 01 

06 300*000 5000.000 

08 0 ,455106746-11-0,436805156- 15 

01-0.163816656*02 0.24210316E-05 

05 0.296399496 01 

06 300.000 5000*000 

06 0,137566466-10 0,142245426-15 
01-0,107788586-02 0,968303786-06 
04 0,493700096 00 
06 300*000 5000*000 

06 0.362018586-10*0*289456276-14 
01-0,187821846-02 0,705545446-05 
04 0,430527786 01 
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Table 6-4. SPECIES WITH THERMODYNAMIC DATA PROVIDED 


ALjS) 

BCL 

BEO(S) 

AL(L) 

BCL+ 

BEO(L) 

AL 

6CLP 

BEO 

AL+ 

BCL2 

BEOH 

ALB02 

6CL2+ 

BEOH+ 

ALCL 

BCL2- 

BE02H2 

ALCL+ 

BCL3 

BE20 

ALCLF 

BP 

BE20F2 

ALCLF2 

BF2 

BE202 

ALCL 2 

BF2+ 

BE303 

ALCL2+ 

•BP2- 

§£404 

ALCL 2- 

BF3 

BR 

ALCL2F 

BH 

BR2CL) 

ALCL3(S) 

BHF2 

BR2 

ALCL3CL) 

BH2 

C(S) 

ALCL 3 

6H3 

C 

ALP 

8N(S) 

C+ 

ALF+ 

BN 

C- 

ALF2 

BO 

CCL 

ALF2+ 

BOCL 

CCL2 

ALF2- 

EOF 

CCL3 

ALF3 (S) 

BOF2 

CCL4 

ALF3 (S) 

B02 

CF 

ALF3 

B02- 

CF2 

ALH 

BS 

CF3 

ALN(S) 

B2 

CF4 

ALN 

B20 

CH 

ALO 

8202 

CH2 

ALO+ 

B203CL) 

CH20 

ALOCL 

8203 

CH3 

ALOF 

B303CL3 

CH4 

ALOH 

B303F3 

CN 

ALOK+ 

BE(S) 

CN+ 

ALOH- 

B£(L) 

CN- 

AL02 

BE 

CN2 

AL02- 

BE+ 

CO 

AL02H 

BEB02 

CCCL 

AL2CL6 

BECL 

CCCL2 

AL2F6 

BECL+ 

CCF 

AL20 

BECLF 

COF2 

AL20+ 

BECL2 (S) 

COS 

AL202 

BECL2 (L) 

C02 

AL202+ 

BECL2 

C02- 

AL203(S) 

BEF 

CP 

AL203(L) 

BEF2 (S) 

CS 

AR 

BEF2 (S) 

CS2 

AR+ 

BEF2 (L) 

C2 

B(S) 

BEF2 

C2- 

B(L) 

BEH 

C2CL2 

B 

BEH+ 

C2F2 

B+ 

BEN 

C2F4 


C2H 

FECL2(S) 

KOH 

C2HF 

FECL2(L) 

. KOH (S) 

C2H2 

FECL2 

KOH(L) 

C2H4 

FECL3(S) 

K2 

C2N 

PECL3U) 

K20(S) 

C2N2 

FECL 3 

LI(S/ 

C20 

FEO(S) 

LI(L) 

C3 

F£0(L) 

LI 

C302 

FEO 

U+ 

C4 

FE02H2(S) 

LICL(S) 

CS 

FE02H2 

LICLd) 

CL 

FE03H3(S) 

LI CL 

CL+ 

FE203(S) 

LIF(S) 

CL- 

F£304(S) 

UFCL) 

CLCN 

H 

LIP 

CLP 

H+ 

LIF2- 

CLF3 

H- 

LIFO 

CL 0 

HALO 

LIH(S) 

CL02 

HBO 

likcl) 

CL2 

HBO+ 

LIH 

CL20 

HB02 

LIN 

CS(S) 

HCL 

LIO 

CS(L) 

HCN 

LIO- 

CS 

HCO 

LIOH(S) 

CS+ 

HCO+ 

LIOHtt.) 

attnt.M 

HOP 

UOK 

CSCLC5) 

HP 

LION 

CSCL(L) 

HNO 

LI2 

CSCL 

H02 

LI2CL2 

CSF(S) 

H2 

LI2F2 

CSFa) 

H20(S) 

LI20(S) 

CSF 

H20(L) 

LI20(L) 

CSO 

H20 

LI20 

CS2 

H202 

LI202 

CS2CL2 

H2S 

LI202H2 

CS2F2 

H3B306 

LI3CL3 

CS2C 

HE 

LI3F3 

£ 

HE+ 

MG(S) 

f 


MG CL) 

P- 

K(L) 

MG 

FCN 

K 

MG+ 

PO 

K+ 

MGCL 

P02 

KCL(S) 

MGCL+ 

F2 

KCLCL) 

MGCLF 

F20 

KCL 

MGCL 2 (S) 

FE(S) 

KF(S) 

MGCL2CL) 

F£(S) 

KF(L) 

MGCL2 

PE(S) 

KF 

MGF 

FE(L) 

KF2- 

MGF2 (S) 

PE 

K2F2 

MGF2(L) 

FECL 

KO 

MGF2 
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fable 6-4. (cont'd) 


MGH 

O 

SIN 

MGN 

o+ 

S10 

MGO(S) 

0- 

S102CP) 

MGOCL) 

OH 

SI 02 (S) 

MGO 

0H+ 

SI 02 (S) 

MGOH 

OH- 

SI02CL) 

MGOH+ 

02 

sro 2 

MG02H2 

02- 

SIS 

N 

P 

S12 

NP 

P(S) 

SI2C 

WF2 

P+ 

SI2N 

N?3 

PCL3 

S13 

NH 

PF3 

XE 

NH2 

PFS 


KH3 

PH 


NO 

PH3 


NO+ 

PN 


NO CL 

PO 


NO? 

PS 


NO?3 

P2 


N02 

P4 


N02- 

S(S) 


N02CL 

SO) 


N02F 

s 


N2 

s+ 


N2C 

SF4 


N2H4 

SF6 


N20 

SH 


N204 

SN 


NA(S) 

SO 


NA(L) 

S0F2 


NA 

S02 


NA+ 

S02F2 


NACL(S) 

S03 


NACL(L) 

S2 


NACL 

S1(S) 


NAF(S) 

SI CL) 


NAF(L) 

St 


NAF 

S1+ 


NAF2- 

stc 


NAH 

S1C2 


NAO 

Sid 


NAO- 

SI CL 2 


NAOH(S) 

std3 


NAOH(L) 

SI CL 4 


NAOH 

sir 


NA2 

SIF2 


NA2CL2 

S1F3 


NA2F2 

SIF4 


NE 

S1H 


NE+ 

S1H4 
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6.1.1 THERMODYNAMIC DATA BELOW 300°K. 

latl^TuTv tT** d ”* *** U " 0t re,ulrM ' However > for low temperature oalou- 

£11. (see se-tlon CUm> fit ^ *» Wernod^c Data 

' e section 6«1). The lover temperature limit t 

'applied with the program la 300°K. *’ Data 

Oiehnodynamlc Data below the temperature T 
ascribed below. temperature, T £ , nay be input by data cards as 


card 1 
card 2 


-ard 3 


L0W T CPHS 


n 


oooo 

T i K < V T ®r 1 


SSTSm.?? •f*J t « f low te«p.ratura 
CRiS tables (col 1 through 10). 

12 character species name, left justi- 
fied, followed by the integer n 
punched in coluan 21. Tiileri 
SJf* such that 1<n<3 ahd represent 
-Ole number of Thennodyanmic Data 
points to be input for this species. 

First Thermodynamic Data point for the 
above species, input 4F 10.0, 15. 


card n+2 


(final 

card) 


^ % nth ( 1<n< 3) nth Thermodynamic Data point for the 

above species, input 4F10.0, 15. 

Repeat cards through n.2 above for each species to be Input. 

Temperature must be T <T <T <T 

12 3 1 * 

END L0W T CPHS end directive (col 1 through 14) 


6-10 




An example of this input is given in Table 6-5 which shows a card listing extending 
the Thermodynamic Data for an 0 2 /H^ propellant to 100°K. Data in Table 6-5 is 
taken directly from the JANAF tables (Reference 23), except for Argon which is 
taken from NASA SP-3001. 

The quantity hJ is defined as 

- (HO,4i 298 )+ ah ? 29Q » cal/mole 
and 

w, , cal/mole - deg K 
*T 

s£ , cal/mole - deg K 


Ref. 23. Stull, D.R. , Prophet, H., et al., JANAF Thennocheraical Tables, Second 
Edition, NSRDS-NBS 37, National Standard Reference Data Series, National 
Bureau of Standards, June 1971. 
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»BLE 6-5. LOW TEMPERATURE , H» r S", DATA POR AN 0,/H., PROPELLANT 


LO* f Cphs 

4R 

100,0 
200,0 
H 

100.0 
200,0 
H2 

100,0 
200,0 
H2o 
100,0 
200.1 
N2 

100. 1 
200,0 
e 

100.1 
200.1 
OH 

100,0 
200,0 
02 

100.0 
200*0 
END lO" 


.0 

>0 


► 0 
„ ft 


4.9681 

4.9681 

2 

-984,5 

-487,7 

2 

31.556 

34,999 

4.968 

4.968 

51118. 

51614, 

2 

21.965 

25.408 

5.393 

6.518 

-1265, 

•462,0 

24.367 

26,520 

7,961 

7,969 

-59376,9 

-58581,9 

36.396 

41.916 

7.074 

6,989 

2 

-1387, 0 
•68* «. 

2 

38*113 

42.986 

S . 666 
5,434 

58479. 

59036, 

32.466 

36,340 

7,567 

7,309 

6.956 

6.961 

C p HS 

2 

7879, 

6623, 

2 

-1381. 

-685. 

35.852 

41.021 

41.395 

46*218 


1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 - 
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6.2 


Title Cards 


This Input permits labeling of runs with alphanumeric information. 4s many 
title cards as desired may be input in sequence. Card format is as follows: 

col 1-5 col 6-77 

TITLE any alphanumeric information 

It is not necessary to input title cards. 
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' 5.3 


DATA Directive and SDATA Name Hat Innut. 


Ito DATA directive and the $DATA fewllmt input aet described below oust always 
be input. It is required lor all problems since it contains the input that cont- 
rsls which calculation modules are to be executed# 

Ihe first input item must be a single card, called the DATA directive card. 

‘ f0raat 0f this cari 18 as Allows; the letters DATA must be punched in colums 

through 4# The DAT, , card is used to ini ora the program that the $DATA namelist 
:n put is to follow. 


The card following the DATA card must contain the name SDATA, and all cards in 

nanelist in P ut ®et **t start in column 2 or greater. Since Namelist input is 

-*rd interpretive, items can be input in any order. The last card in the set must 
contain SEND. 


Users unfamiliar with Namelist input are referred to their FORTRAN reference 
manual. 


Description Units 

DATA directive card 

Namelist name, read in Subroutine 
PR0BLK 


-•'em 
UTA 
DATA 





5.3.1 Specification of Modulea to be Executed . 

If a module is to be executed, it ia necessary to Indicate the fact by input of 
a module flag as described below. For example, if a problem requires that the 0DE 
module be run, it is necessary to input 0DB ■ 1. Chly certain combinations of 
modules are allowed. These are described in Table 6-6. The module flags are: 

Assumed 


Item 


Description 

Units 

Value(s) 

.•DE 

M 

Set 0DE ■ 1 if the ODE module is to be 
executed. 

none 

0. 

0DK 

a 

Set 0DK > 1 if the ODK module is to be 
executed. 

none 

0. 

IDE 

a 

Set TDE a 1 if the TDE option of the 
M0C module is to be executed 

none 

0. 

TDK 

a 

Set TDK a 1 if the TDK option of the 
M0C module is to be executed 

none 

0. 

TDF 

a 

Set TDF a 1 if the IDF option cV the 
M0C is to be executed 

none 

0. 

TDKIL 

a 

Set TDKIL » 1 if the TDK option of the 
M0C is to be executed with an initial 
data line input as described in 
Section 6.9>2. 

none 

0. 

3LM 

a 

Set BLU a 1 if the BLM option of the 
M0C module is to be executed 

none 

0. 

PFO0PT 

a 

Set PPG0PT • 1 if the constant proper- 
ties option cf the MX module is to be 
executed. See Section 6.9.1 for input 
instructions. 

none 

0. 


Table 6-6: Usage of the Module Flags 


0DE-1 , 


0DE-1, 0DK-1, 
SDK-1 , 


DDE-1, 0DK-1, TDK-1, 


0DE-1, 0DK-1, TDK-1, BLM-1, 

0DE-1, 0DK-1 , TDK-1, BLM-1, 
IRPEAT-" , 

0DB-1, 0DK-1, TDK-1, BLM-1, 
IRPEAT-2, 

IDKIL-1 , 

7DE-1, 


“DE-1, BLM-1, 

IDE-1, BLM-1, IRPSAT-1 or 2, 


Mode of Execution 

0DE is huh alone. Options other than the 
rocket (RKT-T) option are allowed. See 
Section 6.4.3 and Reference 9. 

0D£ is run. 0DK is run with 0DE providing 
start conditions. See Section 6.5.1. 

0DK is run alone With initial conditions 
supplied by the user. See Sections 6.5.1 and 
6.5. 1.1. 

0DE is run. 0DK is run with 0DE providing 
start conditions. TDK is run with 0DK and 
TRANS providing M0C initial data line 
conditions. The number of 0DE and 0DK runs 
will be equal to NZ0NES. 

As above with a BLM run added. 

As above with 0DE, 0DK, and TDK repeated 
using the BLM results. See Section 6.8.4. 
This option 's for including the BLM results 
with regen cooling. 

As above except that the 0DE and 0DK runs a~e 
not repeated because there is no heat input 
correction due to regen cooling. 

TDK is run with a M0C initial data line 
input as described in Section 6.9.2. 

0DE will be run for N20NES. TIE is run with 
wJE chemical equilibrium gas properties, and 
with TRANS providing M0C initial data line 
conditions. — 

As above with BLM. 

As above with TEE repeated. 
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Table 6-6 : Usage of the Module Flags (continued) 


Input 

0DE*1, 0DK*1 , TDF=1 , 


0DE*1 * 0DKs1 , TDF*1 , BLM=1 , 

0DF=1 , ODKsI , TDF=1 , BLM=1 , 
IRr'EATsI or 2, 

BLM=1 , 


PFGOPTsI , 


Mode of Execution 


0DE is run. 0DK is run with 0DE providing 
start conditions. However* the 0DK run will 
have a frozen chemical composition. TDK will 
also be run with a frozen composition. 

As above with BLM. 

As above with TDF repeated. 


BLM is run alone. This option requires a 
large amount of input that is provided 
automatically when BLM is run after TDE, TDK 
or TDF. 

The TRANS and M0C modules are run with 
constant gas properties. See Section 6.9.1* 


6.3.2 Inputs for Control for the Program . 


Assumed 

Item Description Units Value(s) 

IRPEAT = Set IRPEAT * 1 or 2 to request that a none 0 

TDK (or TDE) calculation be automati- 
cally repeated after the BLM module 
has been used to calculate a displaced 
nozzle wall. 

If IRPEAT = 1 , the ODE and ODK module 
executions will be repeated with 
adjusted enthalpies for regen cooling, 
see Section 6.8.H. 

If IRPEAT * 2, the ODE and ODK module 
executions will not be repeated. The 
system enthalpy will be unchanged. 

IRSTRT = The program allows for a limited 
restart capability. If the MOC module 
has been run successfully and units 11,1*5,(23 if TDE), 
and 29 have been saved, then the run 
can be continued by input ing 

IRSTRT = 2, 

This will restart the program after 
the MOC calculation and before the BLM 
calculation. The restart handles all 
options involving BLM. 
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i.% * J 


Description 


Number of zones, N, to be used in the 
anlysis. The ODE and-ODK modules will 
be executed N times. Zone 1 repre- 
sents the flow adjacent to the nozzle 
axis, and Zone N represents flow adja- 
cent to the nozzle wall. 

If SI * 0, Eftgliih units are to be 
used for input and output. If SI * 1 , 
SI units are to be used for input and 
output. The SI units required for 
input are shown in parenthesis. 

_When BLM is run with the MOC module, 
"values defining the boundary layer 
edge conditions; x e , y g , U e , T g , and 

P are automatically calculated and 

stored in the XIN£), RIN0, UE 0, TE0, 
and PE0 arrays. If IOFF is not input, 
then these values are stored starting 
with the first entry in each array. 

If I0FF is input then these values are 
stored starting in the I0FF ♦ 1 entry 
of each array. Values must then be 
input into entries 1 through I0FF. 
The boundary layer calculations will 
start at position XIN0( > ) of these 
arrays using the $BLM input. This 
input allows the user to account for 
the development of the boundary layer- 
in the nozzle chamber, upstream of 
ECRAT. I0FF <100. 

If I0FF is input, do not input NXIN0 
in BLM* 


O' 


. 3,3 Specification of nozzle geometry. 

TO eliminate redundancy, all geometry Inputs tbat are common throughou- the 
nodules are Input here. Dote that «* center of- the nozzle geotet^ coordinate 
ayaten, is at the centerline at the throat plane, and that all coordinates are 
normalized by the nozzle throat radius, r t . Hence, axial positions upstream 0 
throat are always negative numbers. In the figures describing the geometry, posi- 
tive angles are shown as counter-clockwise, and negative angles are shown as clock- 

wise. 


Geometric area ratios at which 


0DE and (DDK print out is to be made are 


specified using the input arrays ASUB(1) and ASUP(I). 


The nozzle geonetry is defined in Figure 6-1. The (SDK calculations 

the downstream end of the combustion chamber with a subsonic area ratio 

as Shown. The circular arcs RI and RWU cannot overlap, thus, it is 
that KRAI, Rl, RWW, ^ 55,8141 68 lnput such that 


start at 
of BCRAT, 
necessary 


/k£aT > 1 ♦ (RI + RWTU) (1 - cos THETAI). 


If this condition is not met, subroutine PRES of OTK will print the terminal error 
message: 


inlet geometry incompatible with initial conditions. 

i, addition, the transonic analysis requires that a value of RWTO > -5 be 

input. 

Tie wall geometry downstream of the nozzle throat can be specif led using any 
one of several options. All of these geometries begin with a circular arc 0 ra- 
lius HWTD which extends from the throat point through an angle of WBTA. Geomet- 
ries that can be Input to both the (JDK and M0C nodules are described in Section 

6.3.3.1* 
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:tera 


Description 


Assumed 
Units Value(s) 


RSI 

m 

Nozzle throat radius, r^ 

in, 

(meters) 

0. 

ASUB(1) 

n 

Subsonic area ratios at which informa- 
tion will be printed. 

none 

0. 

NASUB 

a 

Number of entries in the ASUB array 

< 50. Entries must be nonotoni.' decreasing 
ii'. value. 

none 

0. 

ASUP(1) 


Supersonic area ratios at which infor- 
mation will be printed. Entries must 
be nonctonic increasing in value 

none 

0~ 

NASUP 

. m 

Number of entries in the ASUP array 
< 50. 

none 

0. 

ECRAT — 

-« — 

— Nozzle inlet contraction ratio for use 
in 0DE and 01 SC calculations. 

none 

0. 

RI 

a 

Normalized inlet wall radius. 

none* 

0. 

IHETAI 

a 

Nozzle inlet angle. 

degrees 

0. 

HWTU 

a 

Upstream normalized wall throat radius 
RWTU_> .5 is required. 

none* 

0. 

I TYPE 

s 

Type of nozzle wall to be input. 

I TYPE » 0, if the real wall contour is 

none 

0. 


Input. 

I1YPE ■ 1 * if the potential flow wall 
contour is input. 

If IftPEAT • 1 or 2, then the nozzle 
wall will be displaced by ±$* as cal- 
culated by BLM when Idle TDK (or TDE) 
calculations are repeated. This dis- 
placement is -j* for the real wall 
(I TYPE ■ 0), or +5* for the potential 
flow wall ( HYPE « 1). 


* Normalized by the throat radius, r t 
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6.3»3»1 Exhaust Nozzle Geometry Specification, for 0DK and the M0C Modules 


Item 


Description 


Assumed 
Units Value(s) 


IWALL 


option flag for specifying the 
downstream wall. 

1 cone option (input RWTD, THETA t and 
EPS.) 

2 parabolic nozzle contour option 
(input RWTD, THETA, RMAX, ZMAX) 

3 circular arc nozzle contour option 
(input as for IWALL=2) 

4 nozzle contour (spline) option (in- 
put RWTD, THETA, THE, RS, ZS* NWS) 

5 cone with specified end point, 
(input RWTD, RMAX, and ZMAX) 

6 skewed parabola option. (input 
RWTD, THETA, RMAX, ZMAX, THE) 


none 0 


The items required for the various IWALL options are: 


- -TT\ 

s 

downstream wall throat radius of cur- 
vature ratio** 

none* 

THETA 

= 

nozzle attachment angle 


THE 

s 

nozzle exit angle (input if IWALL*4, 
or 6) 

deg. 

EPS 

m 

nozzle expansion ratio (input if 
IWALL* 1 only) 

deg. 

TMAX 

m 

normalized radius at the nozzle exit 
plane (input if IWALL=2 or 3) 

none* 

ZMAX- — •— 

— * — - 

normalized axial position at the 
nozzle exit plane (input if IWALL ■ 
2 or 3) 

none* 


i 


* Normalized by the throat radius r^ 

** If a comer expansion (i.e. Prandtl-Meyer fan) is desired, a value of RWTD • 
.05 is recommended. Experience has shown that values smaller than this £ive tne 
same result but are computationally less efficient. 


Item 


Description 


Units 


Assumed 

Value(s) 


RS(2) s 

table of normalized wall radii down- 
stream of the nozzle tangency point 
(input if IWALLaM). The input tables 
RS and ZS Start with the second entry 
because the first entry is calculated 
automatically by the program, i«e.» 
RS(1)*r T and ZS(1)=z T in Figure 6-1. 

none* 


The wall angle at this position is 
also calculated so that the spline 
contour will be properly joined to the 
nozzle throat contour. 


ZS(2) 

table of normalized axial position 
downstream Of the nozzle, tangency 
point (input if IWALLsM) . 

none 4 

NWS = 

total number of entries in the RS, ZS 
tables. Includes the -irst entry NWS 
£50 (input if IWALLsM) . 

none 

RZN0RM = 

Optional normalizing factor for the 
RS, ZS, RIM, ZIN tables. For example, 
if RS, ZS, and RIN, ZIN were input as 
dimensional numbers, RZN0RM would be 
the throat radius in those units. 

none 

♦ See the inpu 

t variable RZNORM. 
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6.4 0DE INPUT DATA (ALL PROBLEMS SPECIFYING ODE) 


the 0DE Input data described here is exactly as defined in NASA SP- 273, 
Reference 9, except namelists input $!NFt2 and $RJCNP have been combined into 
a single list named $0DE. Any type of equilibrium calculation available with the 
computer program described lh Reference 9 can thus be computed using the $QfDE 
Input data*. In this document, however, only the ftKT option of namelist is 
described. The RKT option differs from that of Reference 9 for problem types 
Other then single zone 0D£. 

the QfDE input data consists of the following input groups: 

1. REACtANTS directive card, followed by up to IS 

data cards, followed by a blank card, 
specifying reactants. 

2. 0MIT and INSERT directives to omit or insert species for 

equilibrium/frozen calculations. 

3. NAMELISTS directive card followed by input 

namelist $0DE specifying input 
case data. 

6.4.1 REACtANTS CARDS 

This set of cards is required for ail WUL problems. The first card in the 
set contains the word REACTANTS punched in card columns 1 to 9. The last card 
In the set is blank. In between the first and last cards may be any number of cards 
up to a maximum of 15, one for each reactant species being considered. The 
cards for each reactant must give the chemical formula and the relative amount 
of the reactant. For some problems , enthalpy values are required* The format 
and contents of the cards are summarized in Table 6-7. A list of some REACTANTS 
cards is given in Table 6- 8 

Relative amounts of reactants. - The relative amounts of reactants may be 
specified in several ways. They may be specified in terms of moles, mole fraction, 
or mole percent (by keypunching M in card column 53) or in terms of weight, weight 
fraction, or weight percent (blank in column 53). 

Relative amounts of total fuel to total oxidants can also be input. For this 
situation, each reactant must be specified as a fuel or an oxidizer by keypunch- 
ing an F or O, respectively, in column 72 of the REACTANTS card. The amounts 

♦These options include TP,HP,SP, TV, UV, or SV problems, Chap.nan-Jouguet de- 
tonation problems, and Incident or reflected shock problems. 
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TABLE 6-7 REACTANTS CARDS 


Order 

Contents 

Format 

Card columns 

First 

REACTANTS 

3A4 

1 to 0 

Any 

One card for each reactant species 
(maximum 15). Each card contains: 




(1) Atomic symbols and formula num- 
bers (maximum 5 sets) 8 

5(A2. F7.5) 

1 to 45 


(2) Relative weight** or number of 
tholes 

F7.5 

46 to 52 


(3) Blank if (2) is relative weight or 
M if (2) is number of moles 

A! 

S3 


(41 Enthalpy or internal energy®, 
cal mole 

F9.5 

54to 62 


(S) Stale: S, L. or G for solid, 
liquid or gas. respectively 

At 

63 


(6) Temperature associated with 
enthalpy in (4) 

F8.S 

64 to 71 


(7) F If fuel or 0 if oxidant 

At 

72 


(8) Density in g cm* (optional) 

f8.S 

73 to 83 

Last 

Blank 




^Program will calculate the enthalpy or internal energy (4) for species in 
the THERMO data at the temperature (6) if zeros are punched >n card 
columns 37 and 38. (See section Reactant enthalpy for additional in* 
fur null ion.) 

^Relative weight of fuel In total fuels or cxldnnt in total oxidants. All 
reactants must be gifen either all in relative weights or all in number 
nf moles. 
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TABU 6-8 LISTING OT SAMPLE REACTANTS CAROS 


REACTANTS 
M 2. 

N ,7806610 .209795AR. 00*662 


100. 0. 0796.15 P 

100. -7,20216*0296,15 0 


REACTANTS 

Ml, H *, cut. 0 4. 

22,06 

-70730, S296.1S 

r 

Cl, H 1.669550 .C312S6S .-008415 

16.56 

•2999.092L296.15 


Alii, 

9,00 

♦0.0 8298,15 

p 

MCI, 0 1. 

,20 

•1*3790, 8296.19 

F 

M2, 01, 

,16 

•66317,4 1296*13 

P 


REACTANTS 
M 2, 

00 . 

100, 

0. 

0296,15 

P 

0 2, 

00 

100.0 

0.0 

6296,15 

0 


REACTANTS 

N 2, H 8, C 2. 

N 2. H *, 

P 2. 


50,0 1273*, 8 L296.15 P ,766 
50,0 12060. 2296,15 P 1.392 
100, -3030,8922 83,2* 0 1,3* 


reactants 

Ul. 

P 2 . 


100, 0. S296.15 P 

100. -3030, 892t 03,2* 0 1,3* 


REACTANTS 
H 2, M *, 

>0. 

12100, 

L298.13 

P 1,093 

Ml, 

70. 

0.0 

8296,15 

P 1.85 

M2, 02. 

170. 

•44080, 

L296.15 

0 1.407 


•Listed above arc six example*. Each example must end with a blank card. 
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given on the REACTANTS cards ore relative to total fuel or total oxidant rather than 
total reactant. 

There ere four options In the $0DE namelist for Indicating relative 
amounts cl total fuel to total oxidant as follows: 

1 . Oxidant to fuel weight ratio (0F is true) 

2. Equivalence ratio (ERATI0 is true) 

3 . Fuel percent by weight (FPCT is true) 

4. Fuel to air or fuel to oxidant weight ratio (FA is true) 

For each option, except 0D£ with NZ0NES»1, the values are given in the 
0FSKED array of $0DE (described in Section 6.4.3). For 0DE with NZ0NES*1, 
the MIX array is used, as described In Reference 9. 
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Reactant enthalpy. Assigned values for the total reactant are calculated 
automatically by the program from the enthalpies of the individual reactants. 
Values for the individual reactants are either keypunched on the REACTANTS 
cards or calculated from the THERM 0 data as follows: 

Enthalpies are taken from the REACTANTS cards unless zeros are punched 
in card columns 37 and 38. For each REACTANTS card with the "00" code, an 
ent .alpy will be calculated for the species from the THERM 0 data for the tempera- 
ture given in card columns 64 to 71. 

When the program is calculating the individual reactant enthalpy for 
values from the THERM0 data, the following two conditions are required: 

1. The reactant must also be one of the species In the set of THERM 0 data. 

For example, NHj(g) is In the set of THERM0 data but NHg(£) is not. Therefore, 

If NH 3 fe) is used as a reactant its enthalpy could be calculated automatically, 
but that of NH 3 U) could not be. 

2. The temperature T must be In the range Tj ow /l . 2 s T s T h h x 1 . 2 where t 

t0 T high 18 temperature range of the THERM 0 data. 9 " ° W 

For cases with N20NES > 1 (see Problem card, Section 6.3) it may be de- 
sirable to modify the enthalpy of each zone. This can be done by using the DELH 
Input array. For the i th zone the I th DELH entry will t* added to the system 
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enthalpy as computed by 0DE from the reactants cards (see above). For example, 
overall system enthalpy of the propellants in the tank can be input through the reactants 
cards and the work adoed or extracted per zone can be input by the DELH entries. 

An. alternate method would be to Input zero enthalpy on the Reactants cards and 
Input enthalpy per zone by the DELH entries. 

6.4.2 0MIT and INSERT Cards 

0MIT and INSERT cards are optional. They contain the names of particular 
species in the library of Thermodynamic Data for the specific purposes discussed 
below. Each card contains the word 0MIT (in card columns 1-4) or INSERT (in 
card columns 1-6) and the names of from 1 to 4 species starting in columns 16, 

31, 46, and 61. The names must be exactly the same as they appear In the 
THERM 0 data. 

6.4.2. 1 0Mr Cards 

These cards list species to be omitted from the THERM 0 data . Tf 0MTT 
cards ar not used, the program will consider as possible species all those 
species .r> the THERM 0 data which are consistent with the chemical system 
being considered. Occasionally it may be desired to specifically omit one or 
more species from considerations as possible species. This may be accomplished 
by means of 0MIT cards. 

6.4. 2.2 IN SERT Cards 

These cards contain the names of condensed species only. They have 
been included as options for two reasons. 

The first and more important reason for including the INSERT card option is 
that, in rare instances, it is impossible to obtain convergence for assigned enthalpy 
problems (HP or RKT) without the use of an INSERT card. This occurs when, by 
considering gases only, the temperature becomes extremely low. In these cases, 
the use of an INSERT card containing the name of the required condensed species 
can eliminate this kind of convergence difficulty. When this difficulty occurs, 
the following message is printed by the program; "LOW TEMPERATURE IMPLIES 
CONDENSED SPECIES SHOULD HAVE BEEN INCLUDED ON AN INSERT CARD". 
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The second and less Important reason is that if one knows that one or 
several particular condensed species will be present among the final equilibrium 
compositions for the first assigned point, then a small amount of computer 
time can be saved by using an INSERT card. Those condensed species whose 
chemical formulas are Included on an INSERT card will be considered by the 
program during tho initial iterations fcr the first assigned point. If the INSERT 
card were not used, only gaseous species would be considered during the initial 
iterations. However, after convergence, the program would automatically Insert 
the appropriate condensed species and reconvene. Therefore, it usually is 
Immaterial whether or not INSERT cards are used. For all other assigned points 
the Inclusion of condensed species is handled automatically by the program. 


6.4.3 jfg-DE NAMELIST INPUT 


The QfDt subprogram contains namelist Input sections $0D£ 
and $SHKINP. the Namelist $0D£ must be preceded by a card with NAMELISTS 
punched in card columns 1-9. 

The S0DE Namelist is required if 0OE “1. or TDK *1, in SDATA as 
descrited in section 6.3. 

For the t-'DE problem type any of nine (9) different equilibrium problems can be 
Solved. These are TP,ftF*SF,lV, W,SV,ftKT,D£TN, and SH0CK. For the 0D£- 
0DK t 0DE-0DK-TDK, or TDE problem, type, only the RKT problem can be 
Solved. In this section only the RKT input option is discussed. Reference 9 is 
•to be used to r repare input for the other equilibrium problems. 

The variables Input by the $0DE namelist are listed in-lable 6-9. 

Additional information about some of these variables follows: 

Pressure units. - Th* program assumes the pressure in the P schedule to be in 
units of atmospheres unless either PSIA * true, or ST = true. 

Relative amounts of fuel(s) and oxldlzer(s) . - These quantities may be specified by 
assigning 1 t: 15 values for either o/f, %F, f/a, or r. If no value is assigned 
for any of these, the program assumes the relative amounts of fuel(s) and oxidizer(s) 
to be those specified on the REACTANTS cards. (See discussion in REACTANTS 
Cards, Section 6.4.1) 

RK T problem . - Only one value for chamber pressure, P, is to be Input for cases 
with NZ0NES > 1 (see $data input. Section 6,3). The stagnation pressure used for 
the l zone will be the value input for P multiplied by the I th value input in the 
schedule XP. -If. not input, all XP entries are assumed equal to one. For TDK type 


problems, zone one is taken about the nozzle axis of symmetry and the lart 
zone is bounded by the nozzle well. Similarly, the I th zone will have a mixture 
ratio equal to the l^ entry In the QfrSKED schedule. 

Print out will be given for the chamber pressure condition (l*e. Stagnation) 
and the throat condition. Print out may be requested at other conditions by use 
of the PCP schedule and the~SUBAR and SUPAR schedules. 

The program will calculate both equilibrium and frozen performance unless 
PR0Z » F or EQL « F are Input. If FR0Z s F, only equilibrium performance will 
be calculated. If EQL « P, only frozen performance will be calculated. 
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T-.3LE 6-9 VARIABLES IN $0D£ NAMELIST 


Variable 

No < of 
entries 

Type 

Value 

before 

read 

Definition and comments 

♦ 

• 

RjCT 

1. 

L 

False 

Rocket problem 3 

P 

26 

ft 

0 

Assigned pressures: stagnation pres- 
sure for rocket problems: values in 
atm unless PSIA, or SI * . T. , (see be- 
low) 

s: 


L 

False 

a Values in P array are in N/m a 

FilA 

1 

L 

False 

a Values in P array are in psla units 

XT’ 

50 

R 

1. 

Multipliers for the i** 1 zone stagna- 
tion pressure (zone 1 * inner zone) 

<? ? . 

1 

L 

False— 

Oxidant to fuel weight ratios are to 
be input 3 

E^A'fl0 

1 

L 

False 

Equivalence ratios are to be input 8 

r?CT- 

1 

L 

False 

Percent fuel by weight are to be input 8 

FA 

1 

L 

False 

Fuel to air weight ratios are to be input 8 

WFSKED 

50 

ft 

0 

For a Rocket problem , and NZ0NES 5 s 1, 
0FSKED will be used rather than MIX 
(see Reference 9). Relative amounts of 
total oxidant to total fuel are input as 
defined by 0F, ERATI0, FPCT, or FA. 
For 0DE-0DK-TDK and TDE problem 
types these values define the oxidant to 
fuel ratios for each zone (zone 1 * 
inner zone) 

DELH 

50 

R 

0 

Corresponding to each zone this value 
will be added to the system enthalpy 
Input thru the reactants cards. Units 
are BTU/# if PSIA=.T. , Joule/kilogram 
If SI=.T . , otherwise Cil/gram. 

:-:lhi 

50 

R 

0 

Corresponding to each zone this value 
will be added to the system enthalpy. 
These values can be used as a 1st est- 
imate for the heat returned to the mair 
combustion chamber by regen cooling 
circuits (0FC input in SBLM). The BLM 
will recalculate these values and, if 
IRPEAT ®1 in SDATA, rerun the ODE, 0DK, 
TDK (or TDE) analysis. Same units as 
DELH. 


' J variable is set to be true. 


Vte: For rocket problems, only one value stagnation pressure can be input. This value 

multiplied by the i th entry of XP will bt used for the stagnation pressure of the 
i th zone. 
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Table f*-y (cont'd) 


Variable 


No. of 
entrie: 


Type 




False 


.0005 


False 


FR0Z 


LISTS? 


False 


KASE 



'if variable Is set to be true. 


Definition and comments 


Consider Ionic species 


Input nozzle mass flow .option for 0DE- 
0DK-TDK or TDt problems. If a value 
for Wfl0W is input an expansion with 
this mass flow will be computed. The 
values input for P and X? are used as 
estimates for computing stagnation pres- 
sure for each zone. The program will 
adjust these stagnation pressures to 
obtain the desired nozzle mass flow 
within a tolerance of RELERR. Units 
are lbs/sec if PSI=.T. otherwise kilo- 
grams/sec. 


Relative difference between requested 
and computed mass flow rate. The pro* 
gram Stops if this error is exceeded * 


Compute and print solutions at these 
values of chamber pressure to pressure 
ratio (entries must be >1.) 


Compute and print solutions at these 
values of subsonic area ratios (en- 
tries must f* 1.1 


Compute and print solutions at these 
values of supersonic area ratio (en- 
tries must / 1) 


Subsonic area ratio to start 0DK cal- 
culations with computed equilibrium 
conditions. The SUBAR input table 
must include an entry equal to ECRAT. 


To start 0DK calculations with com- 
puted equilibrium conditions at the 
nozzle throat. 3 


Calculate rocket performance assum- 
ing equilibrium composition during 
expansion 1 * » 


Calculate rocket performance assum- 
ing frozen composition during expan- 
sion 1 *. 


List names and dates of all species 


residing on thermodynamic data used 


Optional assigned number associated 
with case. 


^Set variable false if these calculations are not desired. 
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6.4.3. 1 Variable Mixture Ratio Option 

The MOC can be run with a variable mixture ratio option b 
retting VARMIXs.TRUE. , arid inputting values into the STREAM ( 1 ) table 
as described below. 

When the variable mixture ratio option is used, there are 
no sliplines in the flow. Instead, the flow mixture ratio will vary 
from the axis (^sO) to the wall ( vsl) as specified in the tables of 
0FSKED vs. STREAM. Sliplines can not be used when the shock option 
is invoked (SH0CKs1 in $M0C), because shock-slipline interaction is 
not provided in the program. Thus, if the shock option is requested 
and there is variation in mixture ratio from streamline to 
streamline, then the variable mixture ratio option must be used. 

The program will not function properly if the spacing in 
the mixture ratio table, 0FSKE3(1), is too large. The required 
spacing depends on the chemical system. As a rule each entry in 
CFSKED(l) must differ no more than 4 or 52 from its adjacent values 
depending on the stoichiometry of the system. There are no spacing 
requirements for the STREAM (1) table. However, the first entry 
must Le 0 and successive entries must increase monotonically wit 
the last entry equal to 1 . 

The tables XP( 1 „ DELH ( 1 ) , and DEEHl(l) of $0DE can be used, 
m which case each entry corresponds to entrie: in 0FSKED(1) and 
STREAM ( 1 ). 

If VARMIXs ..TRUE. , the XM(1) table of *3PK and STRAWS will 
will be computed by the ; rograra and need not be input. 
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Assumed 
Vn) ue(s) 


Item 

VARMIXs 


STREAM ( 1 )s 




De s cription 

If VARMIXs. FALSE. , do not use F 

variable mixture ratio option. 

If VARMIXs. TRUE. , the variable 
mixture ratio option is to be 
used. The following must be 


input ; 



NZ0NES 

in 

$DATA , 

0FSKEDO) 

in 

$0DE and 

STREAM ( 1 ) 

in 

$0DE, below 


If VARMIXs. TRUE., values must 
be input here corresponding to 


the 

0/F 

values input in 

0FSKED ( 1 ) . 


The 

values 

input represent 


(0/F) 

. vs is 1 , . . . , NZ0NES 

where i 

represents the mass 

flow 

between streamline i and 

the 

axis , 

divided by the total 

mass 

flow. 


Thus 

i’ * o 

at the a x ; . , and 


¥ s 1 

at tne wall. 

The 

ith 

entry of 0FSKED is (0/F) i . 

The 

1th 

entry of STREAM is ^ 


0 * 

il i s , . . . , NZ0NES s 50 

and 

V*i.i 
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ff.4.3.2 OPTION TO PUNCH TABLES FOR BOUNDARY LAYER PROGRAM INPUT 
(DOES NOT APPLY TO BLM) 

Conditions computed along the no22le wall can be output as punched cards 
for Input to the BLIMP* TBL, of MABL boundary layer analysis computer programs. 
These conditions are taken by the boundary layer computer program as being the 
lnviscid flow condition at the edge of the boundary layer. Tables to oe punched 
ere: x, y, and P/? c (i.e. the nozzle wall coordinates and the ratio of pressure 
to chamber pressure along the wall). The tables are punched in NAMELIST for- 
mat readable by BLIMP (see the punched card listing given at the end of the sample 
output. Section 7). 

A maximum of SO entries upstream of the throat are saved and punched* 

The wall point at the end of every characteristic is punched up to a maximum of 
500 total table entries • The user may specify a number by which the punched table 
will be offset. Thus, the first point may be output with identification 5 by input 
of 10FF*4. The use of I0FF enables the user to extend the table by adding points 
upstream. 


If punched cards for input to a boundary layer program are required . the 
following items must bo input as part of the $0DE Namelist input: 


Item Name 


Description 


IPTAB 


10FF 

IPtfc, 


If IPTA6=1 . one title card will be punched (this will 
be the last title card input as described in 6.2) 
followed by tables of X, Y. and P/P along the nozzle 
wall* These cards are for Input to vie BLIMP or 
MABL computer programs. The first point punched will 
correspond to beginning of the converging section of 
the nozzle (l.e. at ECRAT; see Figure 6-1. also 
table 6-8). 

The first point to be punched will be numbered as 
10FF+1. 

If 1PTBL end IPTAB-1* tables of M. T/T , C . V. and 
p* Will also be punched* These additional cards are 
required for input to the TBL computer program (i.e. 
the December 21. 1973 version). 


If a TDE problem is specified, the following items must also be input 
When IPTAB=i: RSI, RWTU, THETAI, and RI in $DATA . 


These items define the nozzle geometry from the combustion chamber through the 
throat as shown in Figure 6-1. For a TDE option it is necessary that 1PTBL»0. 
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6.5 jM INPUT DA TA 


0DK input data is required if 0DK *1 or TDK =1 in SDATA as described 
in Section 6.3. The ODK input data consists of three data groups as follows: 


SPECIES 
REACTIONS 
S 0DK 


data group 
data group 
data group 


These data groups are described below in sections 6.5.1, 6.5.2, and 6.5.3, 
resoecti vly. 


6.5.1 SPECIES 

Specie? used by the computer program are determined in several possible 
ways, depending upon the problem type. Methods used to determine chemical 
species for each pfobelm type are discussed below. 

0DK 

For 0DK problems species names and concentrations must be input, see 
Section 6.5. 1.1. 

0DE-0DK 

For CfDE-QTDK problems the initial start conditions for the kinetic expan- 
sion are obtained from an equilibrium calculation. The species list generated by 
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the equilibrium calculation generally contains many more species than the 40 
species for which the ODK subprogram is dimensioned. Therefore a selection 
processes ls-required to interface the ODE calculated equilibrium start conditions 
with the ODK kinetic expansion calculations. This selection is performed using 
the following rules: 

Rule 1 If a species appears In a reaction, it is selected for the kinetic 
Calculation. 

Rule 2 If a species is specified using INERTS directive It is selected 
for the kinetic calculations. 

Rule 3 If any species has a mole fraction greater than an input criterion, 
it is Selected for the kinetic calculation. 

Species which are selected but which dj not appear in a reaction are treated as 

Inert and listed as such on the output list of selected species* 

0DE-0DK-TDK 

For ODE-ODK-TDK problems species are selected by the above rules for 
ODE-ODK problems. However, for multizone TDK cases it is necessary that 
each zone have the same species list. Thus the INERTS input (see Section 6. 5. 2. 7) 
must be used to assure the same species are selected for each zone. 


TDK with Input Initial Line Option 

Species names and concentrations must be input as described in 
Section 6.8.2 when TDKIL -1 is input in $DATA. 

6. 5*1.1 O DK OPTION FOR INPUT OF INITIAL SPECIES CONCENTRATIONS (APPLIES 
ONLY TO THE 0DK PROE1EM TYPEr ~^ 

This Input begins with a single card with SPECIES In columns 1 through 7 
and with either MASS FRACTIONS or MOLE FPuACTIONS in columns 9 through 22. 

If the identifier for mass or mole fractions is omitted, mass fractions are assumed. 
Up to 40 species cards may be input. Only those species specified bv input 
species cards will be considered for an ODK problem . The order of the input 
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species cards is independent of the order In which the species appear on the 
master Thermodynamic Data file. 


A chemical species is identified symbolically by 12 alphanumeric char- 
acters and must correspond, identically with the species name as it appears 
on the Thermodynamic Data file. A complete list of the current species names 
are listed in Table 5-4 (condensed species, however, may not be specified in 
the species list.) The species symbol may not contain the characters * or = . 


Col 

1-10 

11-22 

23-30 

31-60 

61-60 


function 
Not used 

Species symbol (left Justified) 

Not used 

Value of initial species concentration (if zero 
must be input as 0.0) free field P or E format 

User Identification if desired 


6.5.2 REACTIONS 

Chemical reactions must be input if 0 dk »1 or TDK »i in $data. 


Up to 50 reactions with an implied third body and a total 150 reactions may 
by input to the program. Only one card per reaction, and only one reaction per 
card is permitted. Cards specifying third body reactions must precede cards 
specifying all other reactions. Species names appearing in the symbolic reaction 
set must correspond Identically with the species names as they appear in the 
master Thermodynamic Data (see Table 6-4). A card listing for a sample reaction 
set is presented in Table 6-10. 

the symbolic reaction set contains directive cards and reaction/data cards 
in groups as outlined below: 

REACTIONS Directive for start of symbolic reaction 

input 

• Reactions with implied third body speices 

END T3R REAX Directive for end of third body reactions 


All other reactions 


LA£T REAX 
INERTS 

THIRD B0DY REAX RATE RAH0S 


Directive for end of reactions 

Specified Inert Species 

Directive for start of third body reaction 
rate ratios 

Third body reaction rate ratios 


LAST CARD 

The content and format of each 


Directive for end of REACTIONS Input 
type of card is defined as follows: 


6.5.2. 1 The symbolic reaction set begins with a card containing the word 
REACTIONS in columns 1 through 9. Other columns on this card can be used 
lor comments • 


6*5. 2.2 Each card defining a reaction is divided into five fields, separated by 
commas. Each field contains: 

a reaction 

A * followed by a value of A 
N« followed by the value of N 
B * followed by the value of B, 


field 1 
field 2 
field 3 
field 4 


field 5 


the activation energy (Kcal/mole) 
available for comments 


rate parameters 
for the reactioi 


The general form of a reaction is: 


Nj ♦Symbolj + N 2 *Symbol 2 + . . . * N # *Symbol a + *Symbol^ + . . . 

where the left hand side represents reactants and the right hand side represents 
products. The reaction can be either endothermic or exothermic. 

The multipliers, N, must be integers and represent stoichiometric co- 
efficients. If no stoichiometric coefficient is given, the value 1 is assumed. The 
dimensioning currently in the program requires that: 

Nj + N 2 + . . . . s 1° 


and 


N a * N b + • . * * * 10 

The chemical species (denoted by the word "symbol" in the above general 
form) can contain up to 12 characters each of which must match a species name con 
tained in the thermodynamic data (see Table 6-2, card 3). 
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Reaction 

NA++CL-«NACL 

B+2+0-2*B0 

BE+2+2*0H-*BE0H0H 



Na + + Cl’* NaCl 
B^ + OS'" m BQf 
Be*" ♦ 20H~ - Be(0H) 2 


The value assigned to A, N, B define the forward (i.e. left to right) reaction rate 
k, as 

k a a . t -N * e -(1000BAT) 

Iti units of cc, °K, mole, sec. 

All three reaction rate parameters must be Input. The numeric value of each para- 

meter may be specified in either I, F, or E format. If E format is used the E must 
appear before the exponent. 


• 6.5. 2. 3 The reactions with ar. implied third body must precede other types of 
reactions, and must be followed by the directive (columns 1 through 12): 

Col 1 A 

END TBR REAX 

all reactions prior to the above directive will have a third body term added to each 
side of the reaction. E.g. 

H2 * H + H, . . . 

END TBR REAX 

defines the chemical reaction 

h 2 + m*h + h + m 

where M is a generalized third body. Specific third body effects may be Included 

by inputing specific third body reaction rate ratios as outlined in 6. 5. 2. 8. Cards 

encountered after the END TBR REAX directive card do not have a third body term 
added. 

All other reactions are input next, exactly as described under 6. 5. 2. 2. 

6.5.2. 4 After the last reaction has been defined, a card with LAST REAX in 
columns 1 through 9 is input. 
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6.5.2. * Reaction rate data for 13 dissociation-recombination (Implied third 
body) reactions and 35 binary exchange reactions are listed in Table 6-10 for 
propulsion systems containing elements C, Cl, F, N end O. These rates ere 
taken from Reference 16 (two additional reactions from Reference 15 ere Included). 
Cards can be abstracted from Table 6-10 for Input to the computer program. For 
the Implied third body reactions, the third body for which the rate applies is In- 
dicated In parenthesis In the comment field (M represents a "generalized" third 
body, see Section 2.2 for further details). 


6. 5*2.6 INERT SPECIES OPTION 

Inert species (l.e. species not appearing in reactions) can be Included In 
the Input by Input of a card with INERTS In columns 1 through 6 followed by a list 
of inert species names. The species names must e-r-ch be followed by a comma 
and each name must written exactly as in the master Thermodynamic Data. The 
last comma must be followed by the word END. See Table 6-11 for an example. 
The species list can continue on to the next card, but a species name can not 
overlap onto the next card. 


6.5.2. 7 THIRD BODY REACTION RATE RATIOS 


th 


As described above in Section 6.5. 2.2 for the J reaction only one reaction 
rate, kj, where 

kj - at- n j .-y» T • 

can be Input* For three body recombination reactions such as 
H + 0H + M a * H 2 0 + M t 

the rate of reaction Is in general different for each species, depending upon 
the efficiency of the species, M l# as a third body collision partner. As discussed 
In Section 2.2 the temperature dependence of a recombination rate is approximately 
Independent of the third body, l.e. for the I th third body end 1 th reaction: 

ki^AijT-M* 1 

The t'llrd body efficiency of the l 1 * 1 species fe r the J** 1 reaction Is then defined as 
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TABLE 6-10 REACTIONS AND RATE DATA FOR C.CL^H.N, AND O SYSTEMS 
(FROM REFERENCE 16) 


urArtf^s CCLFhno 

i* * a h?a 

ft • w • OH 

0 ♦ * * 0 ? 
r ♦ r • r> 

W * F • Hf 
H ♦ H • *7 
rn ? « o • 

♦ ft i CO 

♦ U ■ w ? 

♦ rt ■ wn 

♦ r • clf 

♦ a ■ hcl 

♦ CL ■ CL* 
EM) OFA* 

M? • ■ 

0M 4 • 

M ♦ f>M B 
ftM B 

Cf> • 

cn • 

Nft • 

#>7 ■ 

QM • 

F? * 

F > 

r? • 
a ? ■ 
o.? « 

M • 

F ■ 

F B 
F? a 
H fl 
H 1 


Cft 


f 

N 

N 

CL 

N 

CL 


ft ♦ 
OH 4 
02 


•MY 


N 
N 
OH 
H 
N? 

H? 

N 
H? 
HfL 4 
HCL 4 
CL? 4 
CL 4 
CLF 4 
CLF 4 


♦ 

« 


clf* h? ■ 

F? • HCL » 

cur * "CL ■ 
f? ♦ a? • 
cn ?* c 


*l.* 197? 
4 * 7 , 5 F 23 
Ah* • OF I* 

Ab|. 2F17 
A-5.7E1S 
4*2, SEl* 
A*6«4F J 7 
4*2. 7F3? 
4*3,«Fl6 
4*1, OEI* 
4*6. AFt* 
A*3.oFl6 
4*3. nEt* 
4*1. IC1* 


JAWNAF 
• N *?.6 
N*1 . 
M»| , 

N*l. 
N*l* 

M4iS 
N*,5 
N*l. 
N*»5 
Na.S 
N*,5 
N*l. 


F9«6 I- j*T**«N 

, N*0.* 

« B*0«* 

• ** 0 «* 

« B*0.* 

• B » 0 «* 

• M 0 i » 

• ■•1*7.5M. 

, B*0«* 

• MO.* 

• 8*0 •* 

• 8 * 0 . * 

* M<)#« 

« Mtn 




h?« 

M?A 

H? 

0? 

cn ? 

cn ? 

n? 

NA 


C ♦ 

c * 
cn ?* 
e • 

HO • 
N • 
Oh ♦ 

H?A* 

hcl • 

PH 


ah 

HO 

N 

n? 

NO 

AH 

F 

r 

AH 

CL 


M • 

0 • 

0 * 

H • 

H « 

0 * 

0 ♦. 
n • 
h?* e? 

HF* F 
HF* H 
HF* HF 
HCu*CL 
HCL *HCL 
h?*CL 
hF*CL 
CfCLF 
i F *CLF 
i hF *CL 
i HA.*F 
hCL*hf 
HF*ClF 
HF.CL2 
CL F *CLF 
Cr>*C^ 
cn * h 
cn* n 
co*no 
Ct* n 

N?*0? 
MA. H 
HF* 0 
hF*oh 
H?0*CL 
HCL*0 


a«2.10C13* 
A>S.75C12* 
ta7.33r.12* 
tal.3Fl3 . 
A * S .* Fl 1 * 
A*«.BSE9 * 
Aa3.}F)3 * 
l«tiOC4 • 
t * l . MEl 3 . 
4*$.3Fl2 • 
AaS.oFi? . 
Aal, 79 £l 0 « 

4 * 3 . OFI * • 

4*1.7NE10* 
As6.2F.11 » 
tal.OFl? . 
4*6. ?Fl? 
As7.6£l2 
tal.BFl? 
Aa5.6El2 
tal.OFlO 
Aal.BElO 
4 * 1 » BE 1 o 
»al»8E} 0 
tal.lEU 
4b5,3F\1 
i taS.SFU 
, 4*1.1211 
, laS«3Fll 
i A*l.0Fl3 
, A.6.3F11 
, ta2.9F)2 
• 4*1 «*Fln 
i A*1.0F11 
. AaS.OEll 


N*0, • Mill'* 

Nan. » M.7P0* 
Nan. * M7,3oO* 
Nao. • 0*0 •* 

Nap. « B*1»0*0* 

F*-.656»MAS.92C» 
Nb£. » Ba» 3 JA. 

N*-l. * 9*6. 2*0. 
Ns.QlS, B*A9«?64» 
Na».$ « BM.OCO* 

Nan. • B*5.7eO* 
NB..S . B«39.739, 
Nan. ♦ 9*3.005* 
n*., 5 • 0*AS.375 
Na»,S « BO.lnO* 

N*»,68« Ba.Oon* 

nbb.of* Ba.Snn* 
N*.«Mi M.300* 
Na.,60* 9*3.200* 
Na-.ns. B*i. 9no* 
Na~»9 « 0atO.337. 
. 0 b 39.A?7. 
. 0aA6.n?%. 
. 8*76,758. 
. 9*6,990* 

. 0*9.6?A* 

. 9»A.3o3» 

. 9*59,61 B, 

. 0*6,55?. 

, b*?9,490. 
. 0*5, 6?**. 
0*.?OO. 


N*-,5 
N*-.5 
N»*,5 
N*«,5 
N*-,5 

N*-»5 

n**,5 

Ne*,5 

N*n. 

N*-.S 


Na»,6*i 
n**. 66, 8*. ton. 
w*»,5 j 0a6,O* 
Ne «.5 ♦ 0 * 0 . 0 * 


0AULCM 


IAR1 NO. 

(AM NO. 

(AM NO. 

(AB) NO. 

(API NO. 

(N?l NC. 

(Ml NO. 

(M NO. 

(N?) NO. 10 
(Ml NO. 11 
(Ml NO. 12 
|M) NO. 13 


1 

2 

3 

A 

5 

6 
T 

• 

0 


BAULCN 

BAULCH 

BAULCN 

BAULCN 

BAULCH 

BAULCN 

BAULCH 

BAULCH 


NO. 21 
NO. 22 
NO. 23 
NO. 25 
NC. A l 
NO. 42 
NO, AB 
NO. 51 
NO. 2* 
NO. 2* 
NO, 27 
NO. 20 
NO. 29 
NO. 30 

no, 2; 

NO. 32 
NO. 2? 
NO, 3* 
NO, 35 
NO. ?• 
NO, 37 
NO. 3B 
NO. 39 
NC. *0 
NO. A3 
NO. ** 
NO, a5 
NO, A6 
NO. a7 
NO. A9 
NC. 50 
NO. 52 
NA, 53 


L»« PF»« 
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thus mjj Is the ratio of the reaction rate with species as the third body to 
the reaction rate input on the reaction card described in Section 6. 5. 2. 2, 

If reaction rate ratios, m^, are to be input for the dissociation-recombin- 
ation reactions, a card with THIRD B0DY REAX RATE RATI0S in columns 1 through 
27 must be input next. If this card ia de leted from the input, the program assumes 
SiL&ijJiLL 2* this card is Included In the input, it must be followed either. by 
a card with ALL EQUAL 1,0 in columns 1 through 13 (which sets all m.,*!) or by 
SPECIES cards as described below: 

The m^, can be input using a card with the word SPECIES in columns 1 through 
7, this word is followed by the name of the i^ species followed by e comma, 
followed by the values in P format, each followed by a comma. These m,. values 
can be continued onto succeeding cards. Note that the m^ values depend on the 
order of input of the reaction cards, i*e. the 1 th reaction is defined by the 
card input after the REAC1X0NS card. 

Table 6-11 gives a sample input for a Hydrogen/Oxygen system using third 
body reaction rate ratios. In this example the three body recombination rates 
are input with Argon as the third body. The rate with H 2 as a third body is a 
factor of 5 larger then with Ar as a third body for the first three reactions -and a 
factor of 4 larger for the fourth (Hydrogen recombination) reaction. 

6,5.2, 9 At this point in the data input deck a card with LAST CARD in columns 
1 through 9 must be input. 
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iSaimnasreaiii 


TABLE «-ll LISTING OF SAMPLE REACTIONS CARDS FOR AN Oj/Hj PROPELLANT 


fcFACtlONS 0*N 

M ♦ OH 0 H?0 
0 ♦ M 0 OH 
0 « 0 * 03 
H * H 0 Hg 
E NO TBA B£A* 


MAY 3-4 1973 
i A*7,5EP3 
* 404,0^10 


4*1 »2El7 

4*6#4£i? 


JAA'NAF 

* N«3,6 
t M»|# 

* n»i. 

* M«i # 


BS#6 

• 0**O.# 

. • *B*0 • t 

• 9®0.f 

i 8*0#« 


H? « Om ■ H # H?0 

OH ♦ OH * 0 ♦ Hpfl 

H ♦ OH 0 0 ♦ HP 

0 ♦ OH * H ♦ Op 

Last R£AX 


• A*3.19E13. N*fl, f B«S.|5* 

• 4«5*7S£iP* N« 9f , 8*.79of 

• 4*7#33ElP* N00# , B*7,3 qC* 

• 4*1 • 3C 1 3 ♦ N« 9# » 0«o.» 


J MEATS N? i AB iENO 

TMtAD BODY REA* RATE RATIOS 

species AP«i..i.«i..!, t 

SPECIES H? » 5 « *5. *S. *4. . 

SPECIES W20 1 30..5.,5, # 3o,, 
SPECIES 03»B.«S,»4*S«i,S# 
SPECIES Mp«4. .4. .4, * 1 ,$« 
SPECIES M»lp«S*i2#S«i2,S*p5*# 
SPECIES 0*l?*5*12*StlP,5*pS.i 
SPECIES OH*i2.5*l2*S»l2,S*35. # 
LAST CARO 


CAR) MO. 
(AB) MO. 
<AP) MO. 
(API MO. 


baulCh mo. ; 

BAULCH NO. J 
BAULCH MO. ; 
BAULCH MO. ; 
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6 . 5.3 S0DK NAMELIST INPUT . 

$0DK Namelist input specifies the conditions for the kinetic expansion calcula- 
tion. The input is read in subroutine 0DKINP and consists of the following groups 
f oata as outlined below: 

6 . 5 . 3.1 Specification of Nozzle Geometry 

6.5*3*2 Integration Control 

6 . 5 . 3«3 Print Control 

6. 5. 3. A Species Selection and Ma3s/Mole Fraction 

Check 


6.5. 3.5 0DK Problem Input 


.5.3.1 Specification of Nozzle Geometry. 


All of the nozzle geometry is to be input using the $DATA Namelist input, see 
Section 6.3*3 and Figure 6-1 . 


An pDK calculation will be carried out for each mixture ratio input in 0FSKED 
:f $0DE whenever 0DK « 1 or TDK ■ 1 in $DATA. 

For a TDK problem, it is necessary that the 0DK calculations be run past the 
nozzle throat. Usually it is not desirable to run the 0DK calculations all the way 
the nozzle exit because of the extra computer time and print out that results 
nowever, if this is desired on a TDK problem, it can be requested by input of item 
EP as described below. 


Description 


Assumed 
Units Value(s) 


Namelist* read in subroutine 
0DKINP. 

If TDK * 1 and a value is input 
here for EP, then eacn ODK will be 
run to expansion ratio EP. 




6. 5 . 3 . 2 INTEGRATION CONTROL 


The integration routine control; tho atop size such that the relative 
error in the dependent variable increments are less than a prescribed fraction, • 
DEL. Only doubling or halving of tho step size is permitted, and on option, 
either all the variables may be considered (JF k Q), or only the fluid dynamic 
variables (J7 a l) may be considered. 


When the flow becomes supersonic and the area defined fluid dynamic 
equations are used, an additional check on continuity is applied in the form 


(pVA)^i - (pVA) N 
(pVA), 


N+l 


< C0NDEL 


where CON DEL is an input relative criterion. 


The step size is held between the two input bounds HMIN and UMAX . 
Fixed step cases may be run by setting input values for HI, HMAX, HMIN all 


equal . 
Item 


Description 

Units 

Assumed 
Va1ue(s ) 

HI 

S 

initial step size 

none 

.01 

HMAX 

S 

upper bound on step size 

none 

0.30001 

HMIN 

s 

lower bound on step size 

none 

.005 

DEL 

= 

fractional incremental 
error 

none 

.001 
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Item 

TEXPLI 

W 

Description 

temperature b- low which 
explicit integration will 
start. Mot recommended, 

Units 

°R 

Assumed 

Value's) 

0 

C0NDEL 

PK 

relative error criterion for 
continuity check for super- 
sonic flow 

none 

C 

O 

* 

IF 

* 0 

all variables considered 
for ’'.top siZ' control 

none 

0 


« 1 

only fluid dynamic variables 

none 



considered for step size 
control, i.c., m,u, and T 


6 . 5 . 3. 3 PRINT CONTROL 

Output from .the Kinetic Expansion Calculation consists of complete 
output for each print station selected. The end point of the nozzle is always 
prir.ed. Print stations are selected from one of the following options: 

Assume' 

Item Description V 3 * ue ^ 


JPRNT 


* -2 print throat and input area ratios (see ARPRNT) -2 

= -1 print at i nternally s ot . a.r_ea_ra_l i_os for conical 

nozzle.* Print at ."ol acted wall contour points 
for contoured nozzles. For the splir.e fit 

option (IWALLsd), prir.t out will occur at 
el:Ch entry in IS of $DATA . For other con- 
tours (IWALLsS or ?. ), print out will 
occur at 20 equally space c axial locations 
along the nozzle. 

* 0 print at every integration step 

a +1 print every ND3rd step up to the throat and 

then nominal area ratios 

s +2 print every ND3rd step over entire nozzle 


*For T PF.NT = -• and a 


ratios arc: 
AkFRNT(l) = 
220, . . .,400 


i / 

* * • t * 


conical .nozzle (i.e. IWAL1 = 1) 

.. , c. ? , 4 0 , • 2 , ■ . . ,58,60, 6*t, ... 
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the internally set area 
16,120,128, ... ,200,210, 


jdSj - 


* 


If JPRNT Is +1 or +2, the following must be input: 


im 

N31 

ND2 

ND3 


Description 

s= first integration step to be selected for 

print 

« last integration step to be selected for 

print 

» print every ND3rd step between ND1 and 

ND2. 


If tprnt is -2, the following must be input: 


Item 


Description 

ARPRNT(I) 

a 

requested area ratios for print, must be 
monotonic increasing and greater than 1.0 
If no values are input, will use 
values from AT"? of $DATA. 

NJPRNT 

s 

number of area ratios requested for print 
< 100. 


An extended print option may be selected as follows: 

Item Value Description 

» 0 no extended print requested 

b 1 extended print option selected (not suggested ) 


DYSCI 


6.5. 3. 4 


SPECIES SELECTION AND MOLE/MASS FRACTION CHECK 


In order to Interface 0DE equilibrium calculated start conditions 
with the kinetic expansion calculations, special consideration must be mode 
for inert species (those not appearing ir the reaction set). Inerts may be select- 
ed explicitly by use of the INERTS directive or by use of a relative selection 
criterion. 

If A MULTIZONE TDK PROBLEM IS SPECIFIED INERTS MUST BE SPECIFIED 
VIA THE INERTS DIRECTIVE. This is required so that the chemistry selected for 
multizone cases will be compatible. 

The INERTS directive is described in Section (6. 5. 2.6). 

The relative selection criterion (0DK or 1 Zone TDK problems,) 

Is described below: 

Item Name 

EPSEL * 


Function 

all species which do not appear explicitly 
ir the reaction set but whose mole fractions 
are greater then the input value for EPSEL, will 
be retained for the kinetic expansion. Species 
selected under this criterion are treated as inert. 
The program assurr.es EPSEL* 1.0E-5, unless input. 


In some instances it may be desirable to use input species concentrations 
which do not sum to unity. Species concentrations, either input or from equili- 
brium start conditions, are summed and the sum checked as described below. 

Function 

Input species concentrations are summed and checked 
versus unity using this input criterion. If 

| 1 ^species concentrations |<XMFTST 

then the test is passed. The species concentrations 
wm then be normalized such that 
^species concentrations * 1 . 

• OE-3, unless input. 

If the test is not passed, an error message will be given and the 
run terminated. 


Item Name 
XMFTST 


The program assumes XMFTST * 
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6.5. 3.5 ODK PROBLEM 1^1 PUT 


This input is required when PROBLEM 0DK is specified on the 
problem card. A kinetic expansion from input arbitrary start conditions is 
to be computed, in addition to the input items described in section 6.5.3, an 
0DK problem requires input of those items described in sections 6.5.1 and 2. 


Item Name 


Iimut Quantity 

Units 

SI Units 

PC 

■1 

chamber pressure 

PSlA 

N/M* 

T 

m 

Initial temperature 

°R 

°K 

V 

m 

initial gas velocity 

ft/sec 

m/sec 

JPFLAG 

•0 

pressure table 
calculated internally 

none 

none 


« i 

pressure table input 


*» 

EC RAT 

* 

initial contraction ratio 

none 

none 

Por JPFLAG * 0 option the following 

must be input: 



Item Name 


Inout Quantity 

UpjtJL 

SI Units 

PI 

B 

initial pressure 

PSXA 

N/M a 

PESTAR 

C 

throat pressure 

PSlA 

N/M 8 

For JPFLAG = 1 option the following must be input: 



Item Name 


In out Quantity 

Units 


PTB(l) 

tr 

normalised pressure 
table entries* 

none 


ZTB(l) 

8 

normalized pressure 
table coordinates** 

none 


NTB 

B 

number of pressure 
table entries, s 127 

none 


2 

6 

initial axial position 

none 



* normalised to input chamber pressure, PC 
** normalized to input throat radius, RSTAR 
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6.S.3.6 MSS AVERAGED 0PK...&E. 

A mass averaged pDK ISP summary page may be obtained at the end of 
the 0DK calculations as described below: 

Item .Name 

MAVXSF "1 Specifies mass averaged ISP option 

XM(1) ■ Ratio of mass flow rate of each zone (zone 1 

■ inner zone) to the total mess flow rate. 


6.6 


$TBANS NAMELIST INPUT. 


When a MPC problem has been apaoified, the Input data set $TRANS is required 
for the transonid calculation. 


Assumed 


Item 


Desarlptlon 

Valuo(s) 

$TRAN5 

■ 

Namelist, read in subroutine TRAN 


XM(1) 

■ 

i 

Ratio of mass flow rate of each zone (zone 1 
■ inner zone) to the total mass flow rate, 
(need not be input if MAVISP • 1 option spe- 
cified and XM input in $0DK). 

50*0 

ALI 

B 

Number of degrees initial line Will be disp- 
laced downstream. The program assumes ALI 
is zero. If ALI is not zero, a symmetric 
throat is required (RWTD s RWTU). 

0 

IBUG 

S 

If input is nonzero, intermediate transonic 
output will be printed. 

0 


The following input may be used to control the construction of the initial line; 


Assumed 

Item Description Value(s) 


MP Number of points to be placed on the initial 50 

line. MP 275. 

A sinusiodal distribution of the following form 
is used: 


r i = [r w sin j)] i = 0,1,2....N 


where N = MP and i3 EXP1 described below. 

Editing is done to control the spacing, see 
DRMIN. 


I ! fill 

EXF1 

DRMIN 


DRMIN1 


Desc ription 


Arsumt ' 

'j-IdaIa) 


* for sinusoidal distribution. 1.2 

p Editing or iter la for oinuaoidal .01 

distribution. The first initial lino point 
below the wall* r 1 , will be spaced auoh that 

r - r. > itin [DRMIN * RWTD, .025] 

W I 


If n v’ilu*. ■?'* :. r ! c .1 .’r ‘o !■.: input in 


a Editing criteria for sinusoidal 5*1 0 

distribution. Points on the initial line 
will be spaced such that 

r - r > DRMIN1 
n n-1 


$END 


6.7 


$MOC NAMELIST INPUT. 


This data set contains the input items for the supersonic Method of Character- 
istics (MOC) module. The items are divided into four types, which are described 
in the following subsections, 

6.7.1 Characteristics Mesh Control 

6.7.2 Inputs from 3DER 

6.7.3 Print Control 

6.7.4 Exit Plano Option 

Often no $MQC input is necessary since the default values are usually suf- 
ficient. 
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a aajti 


6,7,1 Characteristics Mesh Control 


Units 


Assumed 

Values 


Item Name 

Input Quantity 

5M0C * 

Namelist name, read by subroutine 
CHAR 

DS - 

Insertions will be made such 
thet successive points $long 
streamlines will not be sep- 
arated by more than DS. 

DWWl ■ 

insertion control parameter 
A0 W described in Section 5.9.2. 

EPW ■ 

the program will Insert such 
that the wall end point is 
located within a toleranee EPW. 

IMAX - 

the maximum number of Iterations 
to be allowed while attempting to 
achieve e relative convergence for 
the flow van bles of 5*10"5, 

IMAXF - 1 

the program will terminate the case 
if a printed point requires maximum 
iterations for convergence. 

or 

1MAXF • 0 

program will continue the case after 
IMAX iterations per point have occu.-rec 

TEXPLI 

input temperature below which explicit 
integration for the species concentra- 
tions will be used. 

ETHI 

<gfof point editing as described in 
Section 5.9.2, CNTRL. 

ES 

€ g for point editing as described In 
Section 5.9.2, CNTRL. 

DTWI - 

^®t4J crltarioh * or Insertion 111 su * 5 " 

routines INPT, DSP? as described 
in Section 5.9.2, CNTRL. 


none 


degrees 

r 

b • 

none 

.0 

none 

10 

none 

0 


°R, # K if SI Units 

degrees 

.25 

none 

.005 

degrees 

2. 
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■iiirrMi' 


(ji 


• 7 .1,1 M'-C Shook Option 


If the MOC shock option in requested (JH0CK-1), then the 
.rph 'onstructior begins at the 1 r ternect i cn of the initial data 
.it. and the nocsl#* wall, RFiC's are constructed , u '• o r i 7 l 1 


If 


so 


ntiocf j 


i I i»r. free:, t* >- v; ill may or* or,” , 

, f.prrted into the flow field at the oross-over point, Next., l.Hf; 
are constructed starting at the wall, ana the repj on up to t he 
r on n- over point is filled in. The DliC oonetruetj on tre i, continuer. 

• 1th the let point on each IRC bring a right running shook point , 
.hen the axis in reached, the ohock in reflected an a left runnlrp 
shook . The pro graft 1 hen revert » to a RRC eonntruet ion nchome ond the 
hook ] ,o traced until it no chon the wall. It it* then re fleet •. d 
•r cm the wall, and is calculated an a right tunning check 
propagating towards the How axis. Only one shock will be traced, 
■■lit multiple reflections are allowed. 


TEM RACE 
JHCCK* 


1SHCK*= 


A- 


ASSUMED 

INPU T QUANTI TY L 

If 5H0CK*O, shocks will be ignored 
alia the MOC flow field will use LRC 0 

construction. 

If SH0CK«1, a shock will be inserted, 
and traced as described above. 

. f ish :K*1 , no special action taken 1 

If ISHCK-c, crossing of RRCs will be 
ignored for all RRC's that originate 
at the wall upstream of axial coordi- 
nate position XA. If the shock is too 
strong, it cannot be Ignored. 

Axial coordinate position on wall dis- 0 
cussed under ISHCK, above. 




s** 




6.7.2 Flint Control 


Item Name 
N1 


N2 

NC 


MASSFL - 


NDS 


Input Quantity 


flow parameters will be printed 
for every Nl" interior point a- 
long characteristics selected 
for print 

every K2 01 characteristic will 
be selected for print 


for tic / 0 species concentra- 
tions (partial densities) will be 
printed with the flow parameters. 
a NC*1 the quantities A,B,y, 

he !L t (STU/Lb- d R) , and 

enthalpy (ftVsec 2 ) will be append- 
ed to the species concentration 
print. 


at the completion of each left running 

Characteristic (LRC) the messfiow is 

integrated. 

If 

MASSFL ® 0 then no mass flow printed 

MASSFL = 1 then total mass flow and the 
number of points on the LRC 
are printed for each LRC 

MASSFL = 1— then mass flow for each point 
along LRC Is printed 

MASSFL « 3 Same as MASSFL « 2 with the 
addition of execution time at 
the end of each LRC 


see Section 5.8.1, CHAR 


:or NDS « 1 Dividing Streamline Points 
will be printed . (Nor inal) 


?® 8 “ 0 DlvIdin g Streamline Points 
will be suppressed. 


Units 

none 


none 

none 


hone 


none 
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Assured 

Values 

' JOO 

1000 

0 


Jl 


6.7.3 Inputs from DER, Reference 10 


Item Name 
0FBAR 

ETABAR 


DRPISF — « 


Input Quantity 


Overall mixture ratio in- 
cluding condensed phases. 
For print out only. 


Overall evaporation efficiency, 
i.e. the ratio of gas How to 
total propellant flow at the 
throat. K*' total = I 
•ETABAR sp sp,gas 


Unite 

none 


none 


Ratio of total condensed 
phase momentum to the mass 
flow at the throat. Not used, 
reserved for future use. 


Ihf sec/1 bm 
(if SI Units then 
N secAg) 
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6.7.M Exit Plane Option 


On option, the TDK method of charaoteri sties calculation will continue the mesh 
construction through the exit plane of the nozzle and print a summary of the exit 
plane properties. 


Assumed 

Item Description Value(s) 

EXITPL = Exit plane calculation requested if set .FALSE. 

.TRUE, not operational for shock option. 

For the case when two TDK runs are to be — 
made, i.e. when 

IRPEAT *1, or f. in $DATA 

then the exit plane will be computed for the 
second TDK run* but not for the first TDK 
run. 

6.7.5 Punch Initial Line 

During any calculation generating an initial line, the initial line may be 
punched in a form suitable for running an input initial line option. The following 
input is required. 


Item 


Description 


Assumed 

Value(s) 


ILPUCH 

= 

Requests Punching of Initial Line if set 
.TRUE. 

.FALSE. 

IP’JNIT 

- 

The Fortran unit number assigned to the 
PUNCH, (e.g., on Univac 1108 INPUNIT s -3, 
on the other machines it may be 7, 8, etc.) 


I* ’PORT/-' IT NOTE 




If IPTAB s 1 option is selected, i.e., the boundary layer edge conditions punched 
for T3L input, the initial line pinched carps will be interspersed with the TBL 
edge conditions punched cards. 
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$END 


6.8 


Boundary Laver Module (BLM ) Input Data, $BLM 


data jet js required if the BLM is to be executed, i.e., if BLM.1 .0 
«as incut in the $DATA namelist. Most of the data required by the BLM is «-»“■- 

ted « n rrJT.LrTn^ « ^toTL, •**»'- 

Slir!Z tTt id read in will override the assumed or communicated value. 


Ihe input data items to the BLM module are as follows. 


Item 

$BLM 

ITVPE 


WD0T2D 


Description 

Namelist naffle r read in Subroutine 
INPUTB 

flag to specify the type of body 
geometry 

1 , for an aXisymmetric nozzle, re- 
quired by TDK 

2, for an axisymmetric external flow 

3 for a two-dimensional external 
flow starting at a stagnation 
point. 


nozzle mass flow. If the MOC 

nodule was not executed, a value 
can be input here so that a boun- 
dary layer ISP decrement can be 
comDuted 


Assumed 
Units Value (s) 


1 


lb/sec 
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6 - 8.1 


BLM Gas Properties 


Gas properties that are required by the BLM are y, C p , and w versus T, and a 

value for the Prandtl number,. P„. 

r 

If tables for C p and y are not input (see CPO, CKD, and TO, below) then the 
program will prepare these tables using the 0DE module. The tables are prepared 
-sing a series of (T,S) equilibrium calculations, where T varies from 600°R to 
7000°R at 2Q0°R increments. The chamber entropy value is used for S. Values at 
100°R are then extrapolated and added to the table. The table is printed with the 
BLM output. 

If constants defining the gas viscosity, v , and if the Prandtl number, P r , are 
not input, then they will be transmitted from the ODE module. 

If the TD2P and 0DE modules have not been run, then the data listed below must 
be ir^ut. 


Item 


Description 

Units 

Assumed 

Value(s) 

CrCO(l) 

B 

table of ratio of specific heats, 
y , versus T. 

none 

101*1 .4 

: XD 

8 

table of specific heat at constant 
pressure, C p , versus T 

ft/sec-R 

101*6006 

■'XD 

B 

temperatures corresponding to the 
entries in CKD and CPO arrays, 
above 

°R 

300, ...,7000 
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Item 


Description 

Units 

Assure d 

Vslu-r .3) 

NTAB 

S 

number of values entered in CKO, 
the CPO, and TO arrays. 3 < NTAB. 
<101. 

- 

" '1 

RMUI 

* 

Reference viscosity, y , where vis- 
cosity is expressed as 

lbm/ft-sec 

, ?r#- r 



V = U 0 (T/T 0 ) W 



TI 

= 

reference temperature, T Q , for vis- 
cosity. See RMUI, above. 

♦ 

°R 

1500. 

zHega 

s 

Viscosity exponant , u> . See RMUI, 
above . 

- 


PR 

s 

Molecular Prandtl number, P p 

- 

'TC 

. ; - 


0 r:su:v 


va i vh t for 5 ] K,ri 


# , 


ivt Tor air. 


6.8.2 Boundary Layer Ease Conditions 


The coordinates for thy boundary layer are specified in the RIN0 versus XINO 
table. Conditions at the inviscid edge are specifed in the UE0, TE0, and PE0 versus 
XIN0 tables. Conditions at the wall are specified in the TQW and CQW versus X7QV 
tables. The program will redistribute the input stations (up to 201 total) in 
order to have 101 x-stations uniformly distributed per segment, except for the 
first five stations which are generated non-uniformly. The input values of 
y.U e ,T e .P e .fw ^ or Sp* 80(3 ( p V) w are interpolated at the new x-station3 and used in 
the boundary -layer calculations. 

Values of RINO Versus XlNO and the conditions at the inviscid edge of the boun- 
dary layer will be automatically transmitted from the MOC nodule if it has been 
run. Otherwise, they must be input here. See the description of I0FF in $DATA. 
Conditions at the wall must always be input here. 


Item 


Description 

Units 

Assumed 
Value (s) 

XIN0O) 

= 

x a axial coordinate 

none 

- 

R1N0(1) 

s 

y_, radial coordinate 

none 

- 

UE0(1) 

= 

l T e , gas velocity at the inviscid 
edge 

rt /sec 

- 

TE0(1) 

= 

T g , gas static temperature at the 
inviscid edge 

°R 

- 

PE0(1 ) 

s 

P . gas pressure at the inviscid 
edge 

psi 


NXIN0 

= 

number of items in the XINO, RINO, 
720, TE0, and PE0 tables. 

3 < NXINO < 201 




> 3 *' 
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Item 

XTQWO) 

IHFLAG 

TQW(1) 


CQW(1) 

NTQW 


Description 

X , axial coordinate for the TQW 
w 

and CQW tables. 

Flag specifying wall boundary con- 
dition input through TQW array. 

IHFLAG * 0, for temperature 
IHFLAG * 1 , for heat flux* 


T w , wall temperature or wall 
heat flux, depending in IHFLACL- 

For an adiabatic wall, set IHFLAG = 
1 and all TQW(1) * 0. 

For a non-adlabatic wall with 
prescribed heat flux, set IHFLAG* 

1, and note that for heat flux from 
the boundary layer tc the wall . "the 
TQW(l) entries will be negative. 


Units 
none _ 


°R 

or 

B TU 

i n? - s e c . 


(PV) w , mass transfer parameter at lbm/ft^-sec 
the wall 

Number of axial stations, X . Each 

of the above tables must have this 
number of entries. 

3 < NTQW < 201 . 


Assarted 

: ooo.. o. t looo. 

1 

201 * 0 . 


201*0 

3 


TQWC^fs^OI 1 *?.^ 6 S6t f ° r ^ adiabati0 wall » IHFLAG * 1* and 


6.8.3 Integration Step Size Control . 


Item Description Units 


NSEGS s Number of Segments, 1 < NSEGS < 10 none 

The boundary layer wiVT be divfaed into 
segments < f equal length unless values are 
input into below. 

NISPS(I) » Number rf Integration Steps per none 

Segment, £ 101 per segment. 

XSEG(1) * Vector containing the axial (x) lo- none 

cations which define the wall seg- 
ments. The vector is always NSEGS 
♦ 1 Values long. Default values 
are* 

XSEG(i) * z c + (z m -z c ) ( i -1 ) /NSEGS 1*1,2,. . .NSEGS+1 

where and z_ are the end of the 
c m 

cylindrical combustion chamber, and 
the end of the nozzle , 

respectively, as shown in Figure 
6-1 on page 6-20. 

If the boundary layer is to be 
extended upstream of z Q , the usual 

procedure is to input 
XSEG(1) a XIN0( 1 ) 

and XIN0O ) through XIN0(IOFF), 
etc . , of $BLM are input by the user 
as described in $DATA under input 
item I0FF. 

NTR a Station at which transition to none 

turbulent flow is allowed. The 
program starts with an assummed 
boundary layer profile, and then 
turns on the eddy viscosity terms 
for turbulent flow at station NTR. 

For a laminar boundary layer, set 
NTR large, i.e. NTR > n NISPS(I). 


Assumed 

Valu-' - ) 


1 


10*101 


see 

descrip' 

tion 
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6.3.4 Regenerative Cooling Heat Transfer . 


When an engine i3 cooled using a regenerative device, propellant (usually fuel) 
.3 routed around the nozzle so that heat is transmitted from the boundary layer to 
-he coolant. This heat is then returned to the combustion chamber in the form of 
increased propellant enthalpy. When the nozzle wall temperatures are assumed known 
IHFLAG • 0 option), the BLM will calculate the heat flux from the boundary layer 

zo the wall, ^(BTU/ft -sec). These values can also be input directly (IHFLAG > 1 
'Dtion). If the coolant circuit extends from position x g to position x Q and opera- 
tes at an efficiency, n , then the propellant enthalpy entering the chamber will be 
increased by 


AH m 1 nj ft^dA (BTU/lbro) 

*T x o 

where 

dkj, is the total engine mass flow rate 

dA is the nozzle surface area differential, 2irrdx 

Using the method outlined above, the BLM will compute increments of propellant 
enthalpy for up to 3 fuel or oxidizer circuits and print out the resultant enthalpy 
increments. These can be added to a later computer run by using the DELHI (1) input 
' rray. If BLM is to be automatically rerun, then the enthalpy increments will 
automatically be stored into DELHI (1) for the second pass through ODE, ODK, and 
TDK. These enthalpy increments can be calculated in two ways. If the enthalpy in- 
crease is distributed equally throughout the chamber, then 

AH^ • AH 


i . e . 


DELHI (1) - AH, AH, .... 


etc. 


The second method Is to assume that a fuel circuit adds enthalpy only to fuel, 
^nd an o; t circuit adds enthalpy only to ox. It follows that 

:or a fuel circuit 
r+1 

‘ H i • <1yT> 

and for an ox circuit 

r+1 r i 

iH l ■ <lyT> ~ iH - 

where is the mixture ratio of zone i and r is the overall chamber mixture ratio. 

For either method, the steady state engine cycle balance can be approximated as 
follows. First, calculate "adjusted tank enthalpies" for the fuel and for the oxi- 
dizer and input these on the reactant cards. These values must oproximate the 
energy conte:.. of the propellant entering the main combustion chamber accounting 
for all energy gains and losses, except heat returned to the main combustion cham- 
:.er by the regen cooling circuit (s') . Estimates for these amounts are to be entered 
.sing the DELHI (1) input array. An estimate of zero is usually satisfactory. Cor- 
rected estimates will be calculated by BLM and stored in DELHI (1) for a second pass 
through TDK (or TDE). A second pass using these values will be executed automati- 
cally if IRPEAT = 1 was input in the $DATA namelist. 




Description 


XC0(1) 

XCE(l) 

ETAC(l) 

0FC(1) 


DISTRB(I) 


the Ith entry Is the starting position for 
the Ith cooling oircuit*. 

tne Ith entry is the ending position for 
the Ith oooling circuit*. 

the Ith entry is the efficiency for the Ith 
cooling oirouit. 

type of coolant for the cooling circuit: 

0FC(I) * 0. if there is no ith circuit 

0FC(X) si. if the Ith circuit is oxidizer 

0FC(I) s 2. if the Ith circuit is fuel 

I <3 

Flag for method of distributing AH incre- 
ments 


DISTRB(I) s o. for equal distribution of 
heat to chamber 

DISTRB(I) s 1. for distribution of ox heat 
to ox in chamber, and fuel 
heat to fuel in chamber. 


The Ith entry is for the Ith circuit. 


* Normalized by the threat radius, r t » 
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Assumed 
Value! s) 

3*0 

2*0 

3*1 

3*0 

3*1. 




8 • 5 BLM Plotted Output . 


The input described below is used to control plotted 
plot optiona available from BLM are: 

1) Momentum thickness., e, vs. axial position. 

2) Displacement thickness, f\ vs. axial position. 

3) Wall Temperature, T w , vs. axial position. 

M) Verity profiles at specified area ratios, 
locations. 


output from the BLM. Th- 


or at specified axial 


5) Temperature profiles at specified area ratios, or 
locations. 

Example plots are presented in Figures 6-2 through 6-5 for 
above. 


at specified axial 
types 1 through 5, 


Item 

IPR0F 


APR0FC1) 


NPR0F 


IDTPLT 

KMTPLT 

KTWPLT 


Description 

If IPR0F a o, then supersonic area 
ratios are input in APR0F. 

If IPR0F a i , then axial loca- 
tions, x/r t , are input in APRCf 


Area ratios (or axial locations, 

see IPR0F) at which velocity ratio 
and temperature ratio profiles will 
be plotted. Two frames per area 
ratio will be plotted: ’J/U vs. 

. edge 

y/y and T/T vs. y/y 

edge edge edge 

Number of area ratios (or axial 

locations) requested in APR0F. 
NPR0F < 20. 

If KDTPLT a 1 , then displacement 
thickness, £*, vs. axial location, 
x, will be plotted. 

If KMTPLT a i , then momentum thick- 

n fff* 6 * vs. axial location, x, 
will be plotted. 

If KTWPLT a 1 , then wall tempera- 
ture, T w , vs. axial location will 

be plotted. 


Assumed 

Valuers! 

0 


20*0 


0 


0 


0 


0 


■SEND 
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>••9 SPECIAL OPTIONS . 

s*9.1 Constant Properties Gas Option . 

•The TDK program contains a useful option by which the real gas chemistry can be 
replaced by constant properties ** mis try. Multiple zones can be calculated. The 
output includes the nozzle dive . nee efficiency, n DIy , (see Section 5, subroutine 
PRINT). 

The constant properties gas option is run by input of PPG0PT * 1 in $DAIA. 
Sample input data for the constant properties gas option is listed in Table -6-12. 

CW.y the $DATA, $TRAN and $M0C data sets are to be input. These data sets are 
input as described in Sections 6.3* 6.6 and 6.7* respectively, with the following 
required additions to the $TRANS Namelist. 


a(1) * Value of specific heat ratio, y » for each zone, inner to outer* the 

number of zones is specified in $DATA. 

PSA « Chamber pressure in lbs/in 2 . (N/m 2 if SI units) 

*W) * (From Table 6-9) All assumed ■ 1, if not input. 


"( 1 ) 


Chamber temperature, °R, for each zone, inner to outer. (°K if SI 
units) 


:*CC(1) . Real gas constant, ft 2 /sec 2 °R, (i.e., 49721 /M^ . g*J*1 .986/M^ for 

each zone, inner to outer. (m 2 /sec 2 °K if SI unts) 

3W(1) - Qas Molecular Weight. If input then RGC(1) need not be input. 
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Table 6-12: Sample Case for the Constant Properties Cas Option. 


■ ij p 

a 


tlfuE SAM0LE CASE ONC 

data 

$0 A 7 A 
PF30Pf*l, 

*Z0nES«1, 

«SI*2, 

RwTd*2, R*TO=,5, 

THEtA«35,*73a, 

l«AiL*u, 

NwSsti , 

RS(21* 1 *Wi«3, 1,26975, 1,« 79 jo, 1*73575. 2, 0a9tt0,2*a5930,r 
3,6622e>,9**«772,5*79i98,6,J2U<Jt, 

2S(2)s ,39575, *53006, *82905, 1 , 1 9<*73, 1 *66923, 2,32795, 

a,68tl7,7,o8599, 10*9601, 13,31 »<*. 
tH£si 1 *5813# 

SEND 

STRAPS 

691*23, 

PSAtlOO# 

76*5500, 

XMA920, 

X M *i , 

ALt*0, 

SEND 

S*3C 

SEND 
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6.10 


INITIAL VALUES FOR THE $0DK. STRAUS. AND STDK INPUTS 

the following defines nominal values to which variables will be set If 
dot input. If a variable is net-listed, no nominal value Is set. Variables are set 
in the subroutine containing the Namelist read.. 

$0DE, set in subroutine 0DES 


DELS (I) 


0 . 

ecrat 

ss 

0. 

EQL 

e 

..TRUE. 

EQTHST 

BS 

.FALSE, 

ERAH0 

m 

.FALSE 

FA 

* 

.FALSE 

FFCT 

m 

.FALSE, 

FR0Z 

ss 

.TRUE. 

I0NS 

ft 

.FALSE, 

KASE 

ft 

0 

LISTS? 

= 

.FALSE, 

0T 


.FALSE, 

0FSKED(I) 

ft 

0 . 

FC) 

= 

0 . 

PCP(I) 

ft 

0 . 

PSIA 

ft 

.FALSE, 

RELERP 

* 

• 000S 

RKT 

ss 

• FALSE, 

SI 

« 

.FALSE, 

StJBAR(I) 

ft 

0 . 

StlPAR(I) 

ft 

0 . 

WFL0W 

ft 

0 . 

XPd) 

ft 

1. 
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$0DK, set. In subroutine 0DKINP 


C0NDEL 

- 1.0E-6 

DEL * 

* .001 

EPS 

* 0. 

EPSEL 

* 1.0E-S 

HI 

H 

0 
• 

1 

KM AX 

* .10001. 

KM IN 

* .005 

IDYSCI 

= 0 

TWALL 

« 1 

If 

« 0 

JPftAG 

* 0 

JPRNT 

* -1 

TEXPLI 

» 0. 

XMC) 

* 0. 

XMFTST 

*= 1.0E-3 

$ TRANS, set In subroutine TRAN 

ALI 

- 0. 

XBtJG 

* 0 

MP 

o 

cr> 

a 

PMCRlT 

* i . 

PM DEG 

* i. 

XM(D 

* 0. 

$TDX t set It. subroutine CHAR 

DRPISP 

* 0. 

DS 

« .15 

DTWl 

- 2. 

DWWI 

* 3. 

EPW 

“ .01 

ES 

* .005 
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$TDX (coat’d) 


ETA BAR 
ETKI 

xmax 

2MA>a 

MASSTL 

NC 

NDS 

m 

N2 


1. 

.25 

10 

C 

1 

0 

1 

1000 

1000 


'v 





